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FOREWORD 


This  report  was  prepared  by  the  Wave  Dynamics  Division  of  the  Hydraulics 
Laboratory  at  the  U.S.  Army  Engineer  Waterways  Experiment  Station  (WES)  as 
one  of  a  series  of  reports  on  the  General  Investigation  of  Tidal  Inlets  (GITI). 
The  GITI  research  program  is  under  the  technical  surveillance  of  the  U.S.  Army 
Coastal  Engineering  Research  Center  (CERC)  and  is  conducted  by  CERC,  WES,  other 
government  agencies,  and  by  private  organizations.  This  report  contains  de¬ 
tailed  results  of  a  project  on  the  classification  of  inlets  on  the  basis  of 
their  geometry. 

The  report  was  prepared  by  C.L.  Vincent,  formerly  with  WES,  now  Chief, 
Oceanography  Branch,  CERC,  and  W. D.  Corson,  under  the  supervision  of  R.W. 
Whalin,  Chief  of  the  Wave  Dynamics  Division,  and  H.B.  Simmons,  Chief  of  the 
Hydraulics  Laboratory.  CERC  technical  direction  was  provided  by  C.  Mason 
under  the  general  supervision  of  R.M.  Sorensen.  Technical  Directors  of  CERC 
and  WES  were  T.  Saville,  Jr.,  and  F.R.  Brown,  respectively. 

Comments  on  this  publication  are  invited. 


Approved  for  publication  in  accordance  with  Public  Law  166,  79th  Congress, 
approved  31  July  1945,  as  supplemented  by  Public  Law  172,  88th  Congress, 
approved  7  November  1963. 


PREFACE 


1.  The  Corps  of  Engineers,  through  its  Civil  Works  program,  has  sponsored, 
over  the  past  23  years,  research  into  the  behavior  and  characteristics  of  tidal 
inlets.  The  Corps'  interest  in  tidal  inlet  research  stems  from  its  responsibil 
ities  for  navigation,  beach  erosion  prevention  and  control,  and  flood  control. 
Tasked  with  the  creation  and  maintenance  of  navigable  U.S.  waterways,  the  Corps 
dredges  millions  of  cubic  yards  of  material  each  year  from  tidal  inlets  that 
connect  the  ocean  with  bays,  estuaries,  and  lagoons.  Design  and  construction 
of  navigation  improvements  to  existing  tidal  inlets  are  an  important  part  of 
the  work  of  many  Corps  offices.  In  some  cases,  design  and  construction  of 

new  inlets  are  required.  Development  of  information  concerning  the  hydraulic 
characteristics  of  inlets  is  important  not  only  for  navigation  and  inlet  sta¬ 
bility,  but  also  because  inlets,  by  allowing  for  the  ingress  of  storm  surges 
and  egress  of  flood  waters,  play  an  important  role  in  the  flushing  of  bays  and 
lagoons. 

2.  A  research  program,  the  General  Investigation  of  Tidal  Inlets  (GITI), 
was  developed  to  provide  quantitative  data  for  use  in  design  of  inlets  and 
inlet  improvements.  It  is  designed  to  meet  the  following  objectives: 

To  determine  the  effects  of  wave  action,  tidal  flow,  and  related 
phenomena  on  inlet  stability  and  on  the  hydraulic,  geometric,  and 
sedimentary  characteristics  of  tidal  inlets;  to  develop  the  knowl¬ 
edge  necessary  to  design  effective  navigation  improvements,  new 
inlets,  and  sand  transfer  systems  at  existing  tidal  inlets;  to 
evaluate  the  water  transfer  and  flushing  capability  of  tidal  in¬ 
lets;  and  to  define  the  processes  controlling  inlet  stability. 

3.  The  GITI  is  divided  into  three  major  study  areas:  (a)  inlet  classifi¬ 
cation,  (b)  inlet  hydraulics,  and  (c)  inlet  dynamics. 

a.  Inlet  Classification.  The  objectives  of  the  inlet  classification 
study  are  to  classify  inlets  according  to  their  geometry,  hydraulics,  and  sta¬ 
bility,  and  to  determine  the  relationships  that  exist  among  the  geometric  and 
dynamic  characteristics  and  the  environmental  factors  that  control  these  char¬ 
acteristics.  The  classification  study  keeps  the  general  investigation  closely 
related  to  real  inlets  and  produces  an  important  inlet  data  base  useful  in 
documenting  the  characteristics  of  inlets. 

b.  Inlet  Hydraulics.  The  objectives  of  the  inlet  hydraulics  study 
are  to  define  tide-generated  flow  regime  and  water  level  fluctuations  in  the 
vicinity  of  coastal  inlets  and  to  develop  techniques  for  predicting  these  phe¬ 
nomena.  The  inlet  hydraulics  study  is  divided  into  three  areas:  (1)  idealized 
inlet  model  study,  (2)  evaluation  of  state-of-the-art  physical  and  numerical 
models,  and  (3)  prototype  inlet  hydraulics. 

(1)  The  Idealized  Inlet  Model.  The  objectives  of  this  model  study 
are  to  determine  the  effect  of  inlet  configurations  and  structures  on  discharge 
head  loss,  and  velocity  distribution  for  a  number  of  realistic  inlet  shapes  and 
tide  conditions.  An  initial  set  of  tests  in  a  trapezoidal  inlet  was  conducted 
between  1967  and  1970.  However,  in  order  that  subsequent  inlet  models  are  more 
representative  of  real  inlets,  a  number  of  "idealized"  models  representing  var¬ 
ious  inlet  morphological  classes  are  being  developed  and  tested.  The  effects 
of  jetties  and  wave  action  on  the  hydraulics  are  included  in  the  study. 
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(2)  Evaluation  of  State-of-the-Art  Modeling  Techniques.  The  ob¬ 
jectives  of  this  part  of  the  inlet  hydraulics  study  are  to  determine  the  use¬ 
fulness  and  reliability  of  existing  physical  and  numerical  modeling  techniques 
in  predicting  the  hydraulic  characteristics  of  inlet-bay  systems,  and  to  deter¬ 
mine  whether  simple  tests,  performed  rapidly  and  economically,  are  useful  in  the 
evaluation  of  proposed  inlet  improvements.  Masonboro  Inlet,  North  Carolina,  was 
selected  as  the  prototype  inlet  which  would  be  used  along  with  hydraulic  and 
numerical  models  in  the  evaluation  of  existing  techniques.  In  September  1969 

a  complete  set  of  hydraulic  and  bathymetric  data  was  collected  at  Masonboro 
Inlet.  Construction  of  the  fixed-bed  physical  model  was  initiated  in  1969, 
and  extensive  tests  have  been  performed  since  then.  In  addition,  three  exist¬ 
ing  numerical  models  were  applied  to  predict  the  inlet's  hydraulics.  Extensive 
field  data  were  collected  at  Masonboro  Inlet  in  August  1974  for  use  in  evalu¬ 
ating  the  capabilities  of  the  physical  and  numerical  models. 

(3)  Prototype  Inlet  Hydraulics.  Field  studies  at  a  number  of  in¬ 
lets  are  providing  information  on  prototype  inlet-bay  tidal  hydraulic  relation¬ 
ships  and  the  effects  of  friction,  waves,  tides,  and  inlet  morphology  on  these 
relationships . 

a.  Inlet  Dynamics.  The  basic  objective  of  the  inlet  dynamics  study 
is  to  investigate  the  interactions  of  tidal  flow,  inlet  configuration,  and 
wave  action  at  tidal  inlets  as  a  guide  to  improvement  of  inlet  channels  and 
nearby  shore  protection  works.  The  study  is  subdivided  into  four  specific 
areas:  (1)  model  materials  evaluation,  (2)  movable-bed  modeling  evaluation, 

(3)  reanalysis  of  a  previous  inlet  model  study,  and  (4)  prototype  inlet  studies. 

(1)  Model  Materials  Evaluation.  This  evaluation  was  initiated  in 
1969  to  provide  data  on  the  response  of  movable-bed  model  materials  to  waves 
and  flow  to  allow  selection  of  the  optimum  bed  materials  for  inlet  models. 

(2)  Movable-Bed  Model  Evaluation.  The  objective  of  this  study  is 
to  evaluate  the  state-of-the-art  of  modeling  techniques,  in  this  case  movable- 
bed  inlet  modeling.  Since,  in  many  cases,  movable-bed  modeling  is  the  only 
tool  available  for  predicting  the  response  of  an  inlet  to  improvements,  the 
capabilities  and  limitations  of  these  models  must  be  established. 

(3)  Reanalysis  of  an  Earlier  Inlet  Model  Study.  In  1957,  a  report 
entitled,  "Preliminary  Report:  Laboratory  Study  of  the  Effect  of  an  Uncontrol¬ 
led  Inlet  on  the  Adjacent  Beaches,"  was  published  by  the  Beach  Erosion  Board 
(now  CERC).  A  reanalysis  of  the  original  data  is  being  performed  to  aid  in 
planning  of  additional  GITI  efforts. 

(4)  Prototype  Dynamics.  Field  and  office  studies  of  a  number  of 
inlets  are  providing  information  on  the  effects  of  physical  forces  and  artifi¬ 
cial  improvements  on  inlet  morphology.  Of  particular  importance  are  studies  to 
define  the  mechanisms  of  natural  sand  bypassing  at  inlets,  the  response  of  inlet 
navigation  channels  to  dredging  and  natural  forces,  and  the  effects  of  inlets  on 
adjacent  beaches. 

4.  This  report  presents  a  classification  of  inlets  based  on  objective  anal¬ 
ysis  of  similarities  between  inlet  geometric  characteristics.  The  report  con¬ 
tains  substantial  amounts  of  inlet  geometric  data  obtained  from  aerial  photos 
and  boat  sheets  which  may  be  applicable  to  site-specific  studies. 
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CONVERSION  FACTORS,  U.S.  CUSTOMARY  TO  METRIC  (SI)  UNITS  OF  MEASUREMENT 


U.S.  customary  units  of  measurement  used  in  this  report  can  be  converted  to 
metric  (SI)  units  as  follows: 


Multiply 

by 

To  obtain 

inches 

25.4 

millimeters 

2.54 

centimeters 

square  inches 

6.452 

square  centimeters 

cubic  inches 

16.39 

cubic  centimeters 

feet 

30.48 

centimeters 

0.3048 

meters 

square  feet 

0.0929 

square  meters 

cubic  feet 

0.0283 

cubic  meters 

yards 

0.9144 

meters 

square  yards 

0.836 

square  meters 

cubic  yards 

0.7646 

cubic  meters 

miles 

1.6093 

kilometers 

square  miles 

259.0 

hectares 

knots 

1.852 

kilometers  per  hour 

acres 

0.4047 

hectares 

foot-pounds 

1.3558 

newton  meters 

millibars 

1.0197  x  10'3 

kilograms  per  square  centimeter 

ounces 

28.35 

grams 

pounds 

453.6 

grams 

0.4536 

kilograms 

ton,  long 

1.0160 

metric  tons 

ton,  short 

0.9072 

metric  tons 

degrees  (angle) 

0.01745 

radians 

Fahrenheit  degrees 

5/9 

Celsius  degrees  or  Kelvins1 

*To  obtain  Celsius  (C)  temperature  readings  from  Fahrenheit  (F)  readings,  use 
formula:  C  =  (5/9)  (F  -32). 

To  obtain  Kelvin  (K)  readings,  use  formula:  K  =  (5/9)  (F  -32)  *  273.15. 


SYMBOLS  AND  DEFINITIONS 


Symbol  Restriction1 
A  3 

k-Hi  s 

\ a 
AED 


H  1 

B  3 

B.b ij  5 


DCC 


DMA 


DMX 

EC1 


4 
3 

5 


EC2 


EDI 


Definition 

exponent  in  general  regression  equation 

matrix  used  in  formation  of  discriminant  functions 
containing  within-cluster  variance  products 

area  of  the  minimum  width  cross  section 

area  of  the  ebb  delta 

generalized  eigenvector  weighting  coefficient 

exponent  in  general  regression  equation 

matrix  used  in  formation  of  discriminant  function 
containing  population  variance  products 

matrix  containing  distance  coefficients 

coefficients  in  a  generalized  regression  equation 

coefficients  in  discriminate  functions 

depth  at  crest  of  outer  bar  in  channel 

mahalanobis  squared  distance 

average  depth  of  minimum  inlet  width  cross  section 
maximum  depth  of  minimum  inlet  width  cross  section 
depth  at  a  mesh  point  on  a  grid  system 
normalized  depth  at  a  mesh  point 

weighting  coefficient  on  first  eigenvector  of  channel 
profile 

weighting  coefficient  on  second  eigenvector  of  channel 
profile 

weighting  coefficient  on  first  eigenvector  of  ebb  delta 
shape 


ED2 


weighting  coefficient  on  second  eigenvector  of  ebb  delta 
shape 


EMI 

EM2 


weighting  coefficient  on  first  eigenvector  of  M1WC  shape 
weighting  coefficient  on  second  eigenvector  of  MIWC  shape 


EM3 


weighting  coefficient  on  third  eigenvector  of  MIWC  shape 


Restriction  of  a  definition  to  a  given  section. 
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SYMBOLS  AND  DEFINITIONS— Continued 


Symbol  Restriction 


0  -0 
F  ,  f .  . 


jSi > smi 


N,n,nj 


Pi jk ’ Pi 


Definition 


arbitrary  eigenvector 

F  ratio  either  univariate  or  multivariate 

statistical  matrix  of  inlet  data 

coefficient  in  generalized  regression  coefficient 

matrix  of  observations  normalized  to  0  mean  and  unit 
standard  deviation 

arbitrary  function 

mean  value  of  parameter 

coefficient  in  a  generalized  regression  equation 
number  of  clusters 
normalized  depth 

subscripts,  individual  use  varies 
exponent  in  a  general  regression  equation 
subscripts,  individual  use  varies 

standard  deviation  for  variable  computed  within  a 
cluster 

subscripts,  individual  use  varies 
channel  length 

components  of  an  eigenvector 
number  of  inlets  or  variables 
number  of  inlets  or  variables 
number  of  inlets  in  an  analysis 
profile  vector 

probability  that  a  given  inlet  from  a  given  cluster 
belongs  i'n  another  cluster;  posterior  probability 
that  an  inlet  belongs  to  a  cluster 

prior  probability  that  an  inlet  belongs  to  a  cluster 
general  subscript 
coefficient  of  discrimination 


SYMBOLS  AND  DEFINITIONS— Continued 


Symbol  Restriction 

De finition 

r 

number  of  inlets  in  a  given  analysis 

sijk>si  5»6 

discriminant  function 

T,t .  . 

total  variance  matrix  (population  wide) 

u 

ratio  of  pooled  to  total  variance 

Vij 

correlation  matrix 

W  2, 3, 4, 6, 7 

minimum  inlet  width 

W*wij  5 

pooled  variance  matrix 

generalized  coordinate  or  parameter  in  eigenvector, 
cluster,  and  discriminant  analyses 

z*,z 

general  variables  in  eigenvector  analyses 

r 

covariance  matrix 

Y 

dependent  variable  in  a  general  regression  equation 

A 

eigenvalue 

V 

degrees  of  freedom 

£ 

vector  of  normalized  depths  in  an  eigenvector  analysis 
of  the  ebb  delta  shape 

standard  deviation 

4> 

variance  included  in  an  eigenvector  analysis 
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THE  GEOMETRY  OF  SELECTED  U.S.  TIDAL  INLETS 

by 

Charles  L.  Vincent  and  William  D.  Corson 
I.  INTRODUCTION 

1.  Historical  Perspective. 

Along  the  coasts  of  the  United  States,  numerous  inlets  exist  through  which 
water,  sediment,  nutrients,  and  pollutants  are  exchanged  between  the  oceans  and 
estuaries  by  flows  that  are  largely  forced  by  tides.  Among  the  approximately 
300  major  inlets,  a  wide  diversity  in  hydraulic  conditions  and  morphology  is 
found.  Inlet  widths  vary  from  a  few  hundred  to  more  than  20,000  feet  and 
average  depths  vary  from  a  few  feet  to  more  than  SO  feet.  Historically,  the 
economic  importance  of  inlets  has  primarily  been  in  their  role  as  waterways  for 
commercial  navigation,  but  their  potential  for  recreational  interests  has  also 
been  recognized.  More  recently,  the  effects  upon  the  ecology  of  estuarine  areas 
caused  by  changes  in  inlet  configuration  due  either  to  nature  or  man  have  re¬ 
ceived  attention.  The  combination  of  the  large  size  of  the  inlet-estuarine 
systems  and  the  complexity  of  the  physical  processes  makes  the  study  of  inlets 
difficult;  as  a  result,  planning  and  management  of  inlets  often  rely  on  a  less 
than  desirable  information  base. 

2.  Inlet  Classification. 

Both  fixed-bed  and  movable-bed  physical  models  of  tidal  inlets  have  been 
used  extensively  since  1930  to  examine  possible  effects  of  proposed  modifica¬ 
tion  of  inlets  through  construction  of  jetties  and  dredging.  More  recently, 
numerical  models  have  been  formulated  to  simulate  inlet  hydrodynamics.  Both 
types  of  models  provide  valuable  insight  into  the  hydrodynamics  of  the  inlet 
system  but  neither  provide  accurate  quantitative  estimates  of  sediment  trans¬ 
port  or  shoal  and  scour  patterns.  These  models  provide  qualitative  indications 
of  the  effect  of  proposed  modifications  which  are  then  extrapolated  through 
knowledge  of  the  behavior  of  the  prototype  or  other  inlets  to  reach  final  con¬ 
clusions  concerning  the  probable  success  of  the  proposed  modifications.  In  the 
future,  advanced  models  and  additional  knowledge  of  the  physical  processes  should 
improve  the  ability  to  project  the  effect  of  modifications  to  inlets,  but  it  is 
still  likely  that  empirical  knowledge  of  inlet  systems  and  their  response  to 
change  will  continue  to  be  an  important  part  of  the  information  base  for  inlet 
planning  and  management  decisions. 

The  inlet  classification  tasks  of  the  General  Investigation  of  Tidal  Inlets 
(G IT  I j  program  were  formulated  to  improve  the  empirical  data  base  on  inlets  and 
inlet  processes.  It  was  recognized  that  although  a  wide  diversity  of  inlets 
exists,  the  same  basic  physical  processes  occur  at  all  inlets  though  in  differ¬ 
ing  magnitudes  and  settings.  Because  of  these  common  processes,  inlets  which 
appear  to  be  similar  should  be  studied  to  explain  reasons  for  the  similarities. 
The  study  has  several  values.  First,  an  extensive  and  consistent  empirical  data 
base  would  be  collected  for  a  large  number  of  inlets.  Second,  similar  inlets 
would  be  classed  together  to  provide  Corps  personnel  performing  specific  site 
studies  a  guide  to  inlets  with  similar  characteristics.  Third,  the  study  should 
better  define  the  relationships  among  many  inlets,  which  should  provoke  better 
explanations  for  the  similarities  and  differences  found. 
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Early  in  the  planning  phase  of  the  inlet  classification  task,  three  aspects 
of  inlets  were  selected  for  examination:  hydraulics,  geometry  (or  morphology), 
and  stability.  The  efforts  were  designed  to  produce  three  independent  classifi¬ 
cations.  In  a  later  phase  of  research,  efforts  will  be  designed  to  interrelate 
the  three  separate  classifications  and  investigate  the  reasons  for  correspond¬ 
ence  between  well-defined  classifications.  This  report  details  efforts  on  the 
geometric  classification.  Results  from  the  hydraulic  classification  studies 
are  given  in  Jarrett  (1976).  Results  from  the  stability  analyses  are  given  in 
Vincent  and  Corson  (in  preparation,  1980). 

It  should  be  noted  that  all  inlets  are  complex  enough  to  be  considered  dis¬ 
similar.  Yet  as  this  report  shows,  when  the  larger  scale  geometric  properties 
of  inlet  systems  are  analyzed,  striking  similarities  are  found  and  inlet  char¬ 
acteristics  are  shown  to  vary  in  consistent  ways. 

3.  Study  Objectives. 

An  analysis  of  the  morphology  of  inlets  requires  the  study  of  inlet  topog¬ 
raphy  and  bathymetry.  Of  particular  interest  are  the  lengths,  depths,  cross 
sections  and  orientation  of  channels,  and  the  area,  height,  and  location  of 
shoals.  Inlet  morphology  results  from  tide-  and  wave-generated  sediment  trans¬ 
port  over  a  wide  range  of  space  and  time  scales.  The  features  adjust  to  the 
present  hydrodynamic  environment,  but  often  relate  to  past  conditions.  A  fur¬ 
ther  complication  is  the  geology  of  a  particular  inlet  which  may  deform  the 
processes  in  a  unique  way.  Modifications  by  man  can  alter  the  morphology  as 
well.  The  primary  objectives  of  this  study  were  (a)  to  isolate  a  set  of  param¬ 
eters  that  can  be  used  to  quantify  inlet  geometry,  (b)  to  analyze  relationships 
between  the  basic  parameters  selected,  and  (c)  to  analyze  the  relationships 
between  inlets  based  upon  the  parameters  selected. 

The  first  step  in  the  research  is  the  selection  of  parameters  that  satis¬ 
factorily  describe  inlet  geometry.  These  parameters  should  be  representative 
of  the  response  of  the  sediment  mass  to  the  hydraulic  processes.  Additionally, 
for  the  parameters  to  be  useful  in  any  classification  process,  all  must  be 
readily  and  consistently  definable. 

After  consistent  geometric  information  has  been  collected  for  a  large  num¬ 
ber  of  inlets  containing  as  diverse  geometries  as  possible,  the  data  are  ex¬ 
amined  for  relationships  among  the  geometric  parameters.  The  tidal  deltas  and 
channels  develop  as  a  result  of  interaction  between  tide-  and  wave-generated 
flow  fields.  In  the  long  term,  the  magnitudes  of  these  processes  vary  within 
rather  narrow  limits  excluding  the  effects  of  unusually  severe  storms.  Intui¬ 
tively,  it  is  expected  that,  undisturbed,  the  geometry  of  the  inlet  adjusts  to 
the  processes.  Likewise,  even  given  the  varying  influences  of  wave-  and  tide¬ 
generated  phenomena,  consistent  variations  in  the  geometric  parameters  are  ex¬ 
pected  to  occur  representing  the  time-space  integration  of  their  effect.  Few 
such  relationships  have  been  found  previously.  The  best  known  empirical  rela¬ 
tionship  for  inlet  characteristics  is  that  of  O'Brien  (1931),  between  inlet 
throat  cross-sectional  area  and  tidal  prism. 

The  final  task  will  be  classification  of  the  inlets  according  to  their 
geometry.  As  discussed  in  greater  detail  later  in  this  report,  classification 
is  the  process  of  ordering  inlets  on  the  basis  of  the  relationships  among  them. 
Based  on  the  parameters  chosen,  the  similarities  between  inlets  will  be  examined 
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and  the  inlets  organized  accordingly.  The  classification  process  is  difficult 
because  the  inlet  geometry  is  mathematically  multidimensional;  however,  recent 
advances  in  numerical  taxonomy  have  provided  objective  means  for  examining  struc¬ 
tural  organizations.  Thus,  the  differences  among  inlets  can  be  objectively 
measured  and  the  significance  of  the  differences  tested.  The  result  will  be 
a  statistical  basis  on  which  inlet  characteristics  can  be  examined.  Hydraulic 
or  geologic  factors  must  eventually  explain  the  reasons  for  the  organization. 

It  should  be  noted  that  some  of  the  analyses  in  this  report  rely  heavily  on 
the  statistical  methods  of  eigenvector  analysis,  cluster  analysis,  and  discrim¬ 
inant  analysis.  The  motivation  for  use  of  these  methods,  which  have  rarely  been 
used  in  engineering  (but  are  widely  applied  in  the  fields  of  geology,  biology, 
psychology,  and  process  control),  can  be  realized  if  the  futility  of  examining 
manually  the  covariant  properties  of  a  moderate  to  large  number  of  variables 
among  50  or  more  inlets  is  recognized.  The  statistical,  methods  are  required 
to  examine  objectively  the  relationships  among  the  variables  and  the  inlets.' 

The  scope  of  this  report  is  only  to  analyze  the  geometry  of  the  inlet  sys¬ 
tems,  not  their  hydraulics  and  stability.  However,  relationships  which  appear 
to  have  relevance  to  these  other  aspects  of  inlet  systems  are  so  noted. 

4 .  Previous  Classifications. 

A  review  of  the  literature  indicates  that  few  efforts  have  been  made  to 
classify  inlets  according  to  geometry.  Foremost  in  this  effort  has  been  the 
work  by  Galvin,  ^t  al.  (1971a,  1971b)  and  Galvin  (1971)  to  collect  dimensional 
information  on  inlet  geometry  and  relate  inlet  characteristics  to  longshore 
transport  distributions.  This  work  resulted  in  the  definition  of  four  types 
of  inlets  based  on  shape  and  offset  of  the  inlet  flanks:  overlapping  offset, 
updrift  offset,  downdrift  offset,  and  negligible  offset. 

The  other  major  inlet  classification  research  has  been  that  of  Caldwell 
(1955)  in  defining  three  classes  of  inlets  based  on  hydraulic  characteristics. 
Class  I  inlets  have  peak  floodflows  occurring  before  high  tide  by  less  than  1 
hour  with  the  tidal  range  of  the  inlet  and  the  ocean  approximately  equal.  In 
Class  II  inlets  the  peak  floodflow  precedes  high  tide  by  2  to  3  hours  and  tidal 
ranges  in  the  ocean  and  inside  the  inlet  are  about  equal.  The  Class  III  inlets 
have  peak  floodflows  less  than  1  hour  before  high  tide  but  have  tidal  ranges 
inside  the  inlet  substantially  less  than  outside. 

Bruun  (1967)  and  Bruun  and  Gerritsen  (1957,  1959)  related  the  stability  of 
inlets  to  the  ratio  of  the  net  longshore  drift  to  the  maximum  inlet  discharge. 
They  further  suggested  classing  inlets  according  to  their  natural  bypassing 
capabilities. 

The  present  report  is  related  to  the  efforts  of  Galvin,  et  al.  (1971a,  1971b) 
and  Galvin  (1971).  Emphasis  is  placed  on  the  morphologic  characteristics  of  the 
inlets  and  a  quantitative  investigation  of  their  covariant  properties. 

II.  SELECTION  OF  INLETS  AND  PARAMETERS 
1 .  Inlets  Selected  for  Study. 

The  67  inlets  studied  in  this  report  are  listed  by  number  and  geographic 
location  in  Table  1.  The  date  of  the  survey  used,  the  condition  of  the  inlet 


Table  1.  Geographic  order  of  inlets  with  survey  dates, 


inlet  conditions,  and  source  of  data. 


|  Inlet 

|  f.o. 

Name  and  Location 

HBQ3M 

■J 

Source* 

km 

m  am 

i 

Moriches  Inlet,  N.Y. 

mi 

CtGS  5322 

2 

fire  Island  Inlet,  N.Y. 

H  - 

CtGS  7800 

3 

Beach  Haven- Lit tie  Egg  Inlet,  N.J. 

N1 

CI.CS  8220 

4 

Brigantine  Inlet,  N.J. 

1954 

N1 

COGS  8221 

S 

Great  fcgg  Inlet,  N.J. 

1962,  1974 

N1 

GIGS  8676;  NOS  12316 

6 

Corson  Inlet,  N.J. 

1975 

N1 

USAE 

7 

Townsend  Inlet,  N.J. 

1937 

NI 

CIGS  6231 

S 

Hereford  Inlet,  N.J. 

1937 

NI 

CtGS  6236 

9 

Chincoteague  Inlet,  Va. 

1962 

Nl 

CtGS  6764 

10 

Me too* in  Inlet,  Va. 

1862 

NI 

CtGS  794 

11 

Nac hap re ague  Inlet,  Va. 

1911 

NI 

CIGS  3304 

!  '• 

Oregon  Inlet,  N.C. 

1937 

NI 

CtGS  6228 

13 

Hatteras  Inlet,  N.C. 

1916,  1935 

NI 

CCGS  3922  and  5814 

14 

Beaufort  Inlet,  N.C. 

1952 

D 

CtGS  7963 

IS 

Carolina  Beach  Inlet,  N.C. 

1967 

NI 

USAE  CB1  67-7 

16 

Lockwoods  Folly  Inlet,  N.C. 

1924 

NI 

CtGS  44S0 

17 

Shallotte  Inlet,  N.C. 

1934 

NI 

CtGS  5657 

18 

Tubbs  Inlet,  N.C. 

1924 

NI 

CtGS  44S0 

19 

Little  River  Inlet,  S.C. 

1934 

Nl 

CtGS  S6S6 

20 

Murrells  Inlet,  S.C. 

1926,  1974 

NI 

CtGS;  USAE 

21 

North  Inlet,  S.C. 

1925 

NI 

CtGS  4521 

22 

South  Santee  River  Inlet,  S.C. 

1925 

NI 

CtGS  4S22 

23 

Price  Inlet,  S.C. 

1963 

NI 

CtGS  8779 

24 

Capers  Inlet,  S.C. 

1963 

NI 

CtGS  8779 

25 

Dewees  Inlet,  S.C. 

1963 

NI 

CtGS  8779 

26 

Lighthouse  Inlet,  S.C. 

1921 

NI 

CtGS  4181 

,  27 

Stono  Inlet,  S.C. 

1965 

NI 

CtGS  8870 

1  2» 

Fripps  Inlet,  S.C. 

1934,  1972 

NI 

CtGS  5717;  NOS  793 

1  29 

Doboy  Inlet,  Ga. 

1919,  1972 

NI 

CtGS  4099;  NOS  574 

1  30 

Nassau  Inlet,  Fla. 

1934,  1971 

Nl 

CCGS  S79B  and  1110 

!  31 

Fort  George  Inlet,  Fie. 

1924 

CtGS  4376 

32 

St.  Augustine  Inlet,  Fla. 

1924 

CtGS  4453 

:  33 

Ponce  de  Leon  Inlet,  Fla. 

1925 

CtGS  4478 

,  34 

Sebastian  Inlet,  Fla. 

1930 

CtGS  5028 

i  35 

Boca  Raton  Inlet,  Fla. 

1929 

NI 

CtGS  S015 

I  5b 

Hillsboro  Inlet,  Fla. 

1929 

NI  _ 

CtGS  S01S 

j  Gulf  coast  | 

1  37 

Big  Marco  Pass,  Fla. 

1891,  1970 

NI 

CtGS  2038;  NOS  1254 

i  38 

Cordon  Past,  Fl». 

1975,  1970 

J 

NOS  11430  and  1254 

1  59 

Redfish  P.s»,  Fla. 

1960 

NI 

CtGS  8598 

40 

Captiva  Pass,  Fla. 

1960 

NI 

CtGS  8555  and  8362 

41 

Boca  Grande  Pass,  Fla. 

1956,  1970 

D 

CtGS  8358;  NOS  1255 

42 

Gasparilla  Pass,  Fla. 

1956 

NI 

CtGS  8193  and  8196 

43 

Stu^>  Pass,  Fla. 

19SS 

NI 

CtGS  8192 

44 

Mi dn i ght  Pass,  Fla. 

1955 

NI 

CtGS  8154 

45 

Big  Sarasota  Pass,  Fla. 

1954 

NI 

CtGS  8098 

46 

Longboat  Pass,  Fla. 

1953 

NI 

CtGS  8035 

47 

Pass  A  Grille,  Fla. 

1926 

NI 

CtGS  4569 

48 

Clearwater  Pass,  Fla. 

1949 

NI 

CtGS  7875 

!  49 

Pensacola  Pass,  Fla. 

1919 

D 

CtGS  4103 

so 

San  Luis  Pass,  Tea. 

1933 

NI 

CtGS  5488 

51 

Pass  Cavallo,  Tea. 

1934 

NI 

CtGS  5864 

Pac 

ific  coast 

52 

Morro  Bay  Inlet,  Calif. 

1938 

Nl 

USAE  A-2S8S  * 

S3 

Bolinas  Inlet,  Calif. 

1929 

NI 

CtGS  4975 

S4 

Drakes  Inlet,  Calif. 

I860 

Nl 

CtGS  720 

SS 

Bodega  Bay  Inlet,  Calif. 

1931 

NI 

CtGS  5162 

56 

Humboldt  Bay  Inlet,  Calif. 

1S59 

NI 

CtGS  S710 

57 

Coos  Bay  Inlet,  Oreg. 

1885 

NI 

USAE  CB-1'18 

5* 

Umpqua  River  Inlet,  Oreg. 

1903 

NI 

USAE  UM- I'll 

59 

S i us I aw  River  Inlet,  Oreg. 

1891 

NI 

USAE  SL- 1-8 

60 

Alsea  Bay  Inlet,  Oreg. 

1914 

NI 

USAE  AL-1-7 

61 

Yaquina  Bay  Inlet,  Ureg. 

1920 

J 

USAE  Y6-1-63 

62 

Sitet:  River  Inlet,  Oreg. 

1931 

NI 

US>Afc  SE- 1-7/1 

63 

Netarts  Bay  Inlet,  Oreg. 

1957 

NI 

CtGS  8372 

b4 

Tillanmook  Bay  Inlet,  Oreg. 

1910 

NI 

USAE  TM- 1-20 

65 

Nch«i  *  Inlet,  Oreg. 

1957 

j 

CtGS  8368 

66 

Willapa  Bay  Inlet,  Wash. 

1935 

NI 

usaj:  r:-4-:-9i 

Gray'  Harbor  Inlet,  Wash. 

1894 

NI 

USAE  Sheet  524 

_ 

13  Dec.  Report 

**h.n  two  dotes  are  shown,  the  second  date  refers  to  the  survey  used  in  ebb  delta  analysis. 
*NI  indicates  inlet  was  not  jetticd  at  survey  date  nnd  dredging  information  was  w known  at 
survey  date;  0  indicates  inlet  may  have  been  dredged  before  survey  date;  J  indicates  inlet 
has  jetties  in  survey  date  shown. 

’UMI  charts  were  obtained  from  the  area  USAE  Districts. 
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(jettied  or  dredged),  and  the  source  of  the  data  are  given.  About  50  percent 
of  the  inlets  are  on  the  Atlantic  coast,  25  percent  on  the  Pacific,  and  25  per¬ 
cent  on  the  Gulf  of  Mexico. 

The  primary  constraints  used  in  selection  of  the  inlets  and  particularly 
surveys  used  were  (a)  potential  Corps  involvement,  (b)  natural  condition  or 
only  minor  modification  by  man,  (c)  excellent  survey  coverage  of  the  ebb  tidal 
delta,  and  (d)  correspondence  to  inlets  under  study  in  the  stability  analyses. 

As  a  result,  a  wide  range  of  survey  dates  was  used.  Care  should  be  used  in 
applying  specific  data  from  this  report  in  other  analyses  because  later  survey 
data  may  be  available.  It  must  be  emphasized  that  when  an  inlet  is  specifi¬ 
cally  discussed  in  this  report  the  reference  is  related  to  the  condition  of  the 
inlet  at  the  survey  date  which  may  be  very  different  than  the  present  condition. 

More  than  500  boat  sheets  for  more  than  100  inlets  were  analyzed  in  prelim¬ 
inary  studies.  The  inlets  chosen  are  expected  to  typify  the  wide  range  and 
diversity  of  inlet  conditions  without  bias  toward  any  of  the  parameters.  In 
this  way  an  appropriate  sample  of  the  statistical  population  of  inlets  was 
chosen.  Because  a  wide  range  in  dates  and  sources  of  surveys  has  been  used, 
it  is  important  to  recognize  that  the  accuracy  of  the  data  does  vary.  It  is 
expected  that  the  errors  are  random  and  not  biased. 

2.  Orientation. 


Throughout  this  report  the  orientation  of  features  and  graphs  of  the  fea¬ 
tures  will  be  specified  using  the1  convention  that  LEFT  signifies  left  of  an 
observer  standing  on  the  mainland  looking  seaward.  Depths  are  negative  below 
mean  low  water  (MLW) . 

3.  Definition  of  Inlet  Geometric  Parameters. 

To  examine  the  variation  of  inlet  geometry  and  to  derive  an  inlet  classifi¬ 
cation,  it  is  essential  to  select  a  consistently  definable  set  of  parameters 
that  describes  the  basic  characteristics  of  inlet  geometry.  More  than  50  param¬ 
eters  have  been  measured  and  examined  at  various  inlets  to  select  parameters  for 
final  study.  Some  parameters  were  readily  discarded  because  of  definitional 
difficulties.  Others  were  discarded  because  they  proved  insignificant  in  pre¬ 
liminary  cluster  analyses.  The  parameters  selected  for  final  analysis  appear 
readily  definable  at  most  inlets  and  in  previous  tests  proved  to  be  important 
descriptors  of  inlets. 

As  noted  later  in  the  examination  of  the  parameters  used,  there  are  no 
descriptors  of  the  flood  tidal  delta  or  the  bay  system.  Approximately  20  vari¬ 
ables  examined  produced  no  consistent  set  of  meaningful  parameters.  This  is  a 
result  of  the  extreme  complexity  of  the  bay  areas,  the  slow  adjustment  of  these 
areas  to  changes  at  the  inlet,  and  the  great  influence  of  geologic  factors.  The 
bay  channels  result  not  only  from  current  inlet  conditions  but  may  be  remnants 
of  older  tidal  deltas  or  may  result  from  the  influence  of  estuarine  processes 
not  directly  related  to  the  inlets.  It  is  recommended  that  the  characteristics 
of  these  systems  as  related  to  the  inlets  be  considered  as  a  separate  estuarine 
classification  system. 

The  sections  of  the  inlet  system  considered  in  this  report  are  those  parts 
located  beyond  what  is  defined  as  the  minimum  inlet  width  cross  section  (MIWC) 
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and  includes  all  of  the  ebb  tidal  delta.  Characteristics  included  are  the 
dimensions  and  shape  of  the  following: 


(a)  The  MIWC  profile. 

(b)  The  main  channel  from  the  MIWC  to  a  point  beyond  the  crest 

of  the  ebb  tidal  delta. 

(c)  The  ebb  tidal  delta. 

These  are  the  areas  that  are  highly  influenced  by  a  combination  of  waves,  wave¬ 
generated  longshore  currents,  and  tidal  currents.  They  also  include  the  areas 
of  the  inlet  that  frequently  pose  the  most  difficulty  in  maintenance  of  naviga¬ 
tion  channels. 

a.  Channel  Cross  Section  at  Minimum  Inlet  Width  (MIWC).  The  characteris¬ 
tics  of  the  narrow  parts  of  the  inlet  gorge  are  to  a  great  degree  determined 
by  the  exchange  of  waters  between  bay  and  ocean.  The  shape  and  size  of  the 
channel  result  from  a  balance  between  tidal  current  transport  of  sediment 
through  the  channel  and  wave-generated  transport  into  the  inlet  from  adja¬ 
cent  beaches.  Here  attention  is  focused  on  the  channel  cross  section  located 
at  the  minimum  inlet  width  (measured  at  MLW) .  For  hydraulic  considerations, 
the  cross  section  with  minimum  area  has  generally  been  used;  it  is  frequently 
located  close  to  the  MIWC.  Another  consideration  in  the  choice  of  the  MIWC  is 
that  the  minimum  inlet  width  can  be  defined  on  photos  which  are  a  major  source 
of  inlet  stability  data.  The  use  of  MIWC  allows  for  correspondence  between  the 
two  data  sets. 

The  minimum  width,  W,  of  an  inlet  will  be  defined  as  the  minimum  distance 
between  the  shorelines  bounding  the  inlet.  This  definition  can  easily  be  applied 
to  cases  similar  to  the  situation  in  Figure  1.  Fortunately,  this  is  the  dominant 
case  encountered.  However,  two  inlets  similar  to  the  one  in  Figure  2  occurred  in 
the  set  of  inlets  studied.  Figure  2  shows  that  a  minimum  width  near  the  inlet 
throat  may  not  allow  the  parameters  defined  for  this  study  to  be  measured  in  a 
manner  consistent  with  the  other  inlets.  An  alternate  inlet  width  (Fig.  2)  was 
used.  This  width  was  determined  to  be  given  by  a  line  with  minimum  length  par¬ 
allel  to  the  shoreline  bend  and  seaward  of  any  channel  bifurcation.  W  is  given 
in  feet. 

After  the  minimum  inlet  width  was  established  and  measured  at  Mi.  ,  a  detailed 
cross  section  was  plotted.  The  depths  were  recorded  as  a  function  of  distance 
from  the  left.  The  average  depth,  DMA,  was  calculated  and  the  maximum  depth, 
DMX,  noted.  Both  DMA  and  DMX  are  negative  quantities.  The  cross-sectional 
area,  Ae,  was  determined  and  found,  as  expected,  to  be  very  close  to  DMA  *  W. 
Hence,  it  is  a  redundant  parameter.  The  analysis  of  the  cross-sectional  shapes 
is  defined  later  in  this  report.  DMA  and  DMX  are  given  in  feet. 

b.  Channel  Profile.  From  the  MIWC,  the  deepest  channel  was  noted  and 
traced  seaward  (Fig.  3).  The  depths  were  recorded  as  a  function  of  distance 
along  the  length  of  the  channel  from  the  MIWC.  The  channel  was  stopped  at  that 
point  where  minimum  depth  across  the  outer  bar,  DCC,  was  located.  The  channel 
length,  L,  was  defined  as  the  distance  measured  along  the  channel  from  the 
MIWC  to  the  crest  of  the  bar  in  the  channel.  The  analysis  of  the  shape  of  this 
profile  is  discussed  later.  A  definitional  sketch  for  DCC  and  L  is  given 

in  Figure  3.  Profiles  of  the  67  inlets  selected  are  plotted  in  Appendix  A. 

L  and  DCC  are  given  in  feet;  DCC  is  a  negative  quantity. 
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Figure  1.  Measurement  of  inlet  minimum  width  (W) . 


Figure  2.  Measurement  of  alternate  inlet  width. 
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Figure  3.  Measurement  of  channel  parameters  DCC  and  L. 

c.  Measurement  of  Ebb  Tidal  Delta  Area.  Ebb  tidal  delta  area  (AED)  is 
defined  as  the  area  seaward  of  the  inlet  bounded  by  the  shoreline,  the  contour 
depth  at  the  crest  of  the  outer  bar  in  the  channel,  DCC,  to  a  point  where  it 
parallels  the  shoreline  and  the  line  of  inlet  minimum  width  (Fig.  4).  Using  a 
Bruning  Areagraph  chart  No.  4849  and  hydrographic  charts,  the  area  the  ebb 
delta  was  first  calculated  in  square  inches  and  then,  using  the  scale  of  the 
hydrographic  chart,  was  recorded  in  square  miles. 

4.  Representation  of  Channel  Cross  Section,  Channel  Profile,  and  Ebb  Delta 
Shapes:  Mathematical  Considerations. 


An  important  characteristic  of  a  channel  cross  section,  profile,  or  shoal 
is  its  shape.  Unless  the  shape  is  simple,  it  is  difficult  to  express  succinctly 
in  mathematical  terms.  Previous  techniques  have  relied  upon  the  definition  of 
some  simple,  but  arbitrary  parameters  that  can  be  demonstrated  by  the  case  of 
a  channel  cross  section.  A  typical  parameter  selected  to  describe  the  shape 
would  be  the  distance  from  one  side  of  the  inlet  to  the  point  of  maximum  depth. 
However,  if  the  cross  section  had  two  major  channels,  two  such  parameters  would 
be  defined.  As  can  be  seen,  the  solution  is  to  define  N  in  such  parameters. 

If  a  cross  section  has  less  than  N  channels,  the  values  for  the  nonexistent 
channels  would  be  set  zero.  Herein  lies  a  mathematical  ambiguity:  a  zero  value 
implies  a  channel  at  the  base  point  for  measurement.  It  is  possible  to  define 
additional  parameters  to  relieve  the  ambiguity,  but  the  problems  tend  to  be 
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Figure  4.  Measurement  of  ebb  delta  area  (AED). 

compounded.  Another  approach  used  is  to  fit  the  curve  or  surface  with  a  mathe¬ 
matical  function  such  as  a  Fourier  series  or  orthogonal  polynomials.  This  is 
satisfactory  only  if  a  few  functions  are  required. 

Recent  research  has  resulted  in  a  simple  method  for  representing  shapes. 

The  analysis  results  in  a  succinct  parameterization  based  on  rigorously  derived 
mathematical  functions  that,  in  the  sense  of  least  squares  optimization,  are 
those  functions  that  best  fit  the  shapes  under  consideration.  The  method  is 
termed  an  eigenvector  analysis  and  is  directly  related  to  an  R-mode  principal 
component  analysis  (Kendall  and  Stuart,  1968).  It  has  been  used  to  examine 
shape  variations  in  meteorologic  parameters  such  as  spatial  pressure,  tempera¬ 
ture,  and  rainfall  fields  (Kutzbach,  1967);  in  Inner  Continental  Shelf  bathym¬ 
etry  (Resio,  et  al.,  1977);  in  profiles  across  barrier  islands  (Vincent,  et  al., 
1976);  in  beach  profiles  (Winant,  Inman,  and  Nordshrom,  1975;  Vincent  and  Resio, 
1976);  and  in  channel  cross  sections  (Vincent,  1976).  An  heuristic  explanation 
of  the  technique  is  given  in  the  present  report  for  a  channel  cross  section. 

More  rigorous  explanations  of  the  mathematics  can  be  found  in  Kendall  and 
Stuart  (1968),  Resio,  et  al .  (1977),  and  Vincent  (1976). 

A  channel  cross  section  can  be  represented  as  an  M  component  vector,  Z*, 

Z*  =  (dx,  d2,  ...,  dM)  (1) 

where  d^  is  the  depth  at  the  i ^  location  along  a  traverse.  What  is  sought 
through  an  eigenvector  analysis  is  a  series  of  M  geometric  shapes,  represented 
as  M  vectors,  e_  ,  e_2,  £3>  . . . ,  e^  each  with  M  conponents  such  that 

— *  =  H-  +  ai«i  +  a2S.2  +  +  (2) 

where  alf  a2,  ...»  a^  are  a  series  of  coefficients  that  weight  the  individual 
shape  vectors  in  reconstructing  the  original  cross-section  shape  2_  ;  is  a 
constant  vector  representing  the  mean  shape. 
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To  apply  an  eigenvector  analysis  it  is  necessary  to  have  a  set  of  N  cross 
sections,  Z *,  Z^,  Z$,  which  are  typical  of  the  variation  expected  in  the 

shape  under  study  and  with  N  >>  M.  Typically,  these  vectors  are  given  as 
deviations  with  respect  to  the  mean  shape,  jj 


that  is. 


N 

I*  Z  U 

i=  1 


(3) 

(4) 


The  M  x  M  covariance  matrix,  T,  for  the  M  components  of  Z  is  constructed 
from  the  N  observations  of  Z_.  The  shapes  desired  are  solutions  e^  to  the 
matrix  equation 


Te  =  Xe  (50 

where  X  is  an  eigenvalue.  There  are  M  eigenvectors  and  M  eigenvalues 

X  if  r  is  of  rank  M.  The  shapes  derived  in  this  manner  have  several  impor¬ 
tant  properties: 

(a)  For  i  ^  j ,  e^  •  ej  =  0;  the  shapes  are  orthogonal. 

(b)  For  i  =  j ,  e^  .  ej  =  1;  the  shapes  are  unit  vectors. 

(c)  The  magnitude  of  the  eigenvalue  X^  associated  with  shape 
e£  is  expressly  the  variance  in  the  set  of  original  profiles 
explained  by  shape  e^. 

(d)  It  can  be  mathematically  demonstrated  that  by  ordering  the 
shapes  by  descending  value  of  the  associated  eigenvalue,  e  must 
be  that  shape  component  that  explains  the  most  variance  in  the  set 
of  cross  sections;  £2  must  be  that  shape  that  explains  the  next 
most  variance  given  that  is  orthogonal  (or  independent)  of  e, ; 
and  so  forth  for  all  others. 

(e)  Given  an  individual  eigenvector,  e_  ,  with  e_  = 

(l-£l»  •••»  components  represent  the  con¬ 
tribution  of  the  eigenvectors  to  variation  at  transect 

location  j.  The  effect  that  e^  has  in  contributing  to  re¬ 
construction  of  the  individual  profiles  can  be  seen  by  plotting 
1- , ,  1.  ,  ...,  1  as  in  Figures  5  and  6.  (An  analogy  can  be 

1  fy  c.  ifl 

drawn  to  a  harmonic  expansion  where 

n 

f(x)  =  T.  sin  (ix)  : 
i=l 

the  corresponds  to  a^  and  sin(ix)  corresponds  to  e^. 

The  values  l^j  correspond  to  the  individual  value  of  sin(ix) 
for  some  given  x-ordinate  x  • .) 

li 
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Figure  S.  Minimum  inlet  width  cross-section  parameters.  Mean  cross  section  and  pointwise 

standard  deviation  are  shown  along  with  the  shapes  of  the  first  three  eigenvectors 


Thus,  through  an  eigenvector  analysis  shape,  functions  are  derived  that  are 
an  optimal  decomposition  of  shape  variations.  Further,  there  is  a  new  represen¬ 
tation  of  the  profile  in  the  M  dimensional  eigenvector  space;  i.e.,  Z.  is 
transformed  to  a  vector 


a.  —  ( a .  ,  a .  >  •  •  •  >  a . , ,) 

— 1  Z2  iM 


(6) 


where  a are  the  weighting  coefficients  displaying  the  part  of  the  decompo¬ 
sition  of  cross-sectional  shape  explained  by  &_■ . 

The  primary  advantage  of  the  eigenvector  analysis  is  now  apparent.  Having 
ranked  X1,  X2 ,  ....  A^  in  descending  order,  it  is  convenient  to  find  if  there 
is  an  index  K  such  that  A^+1  +  ...  +  A^  is  sufficiently  small  to  be  neglected. 
If  so,  an  estimate  of  the  original  shapes  is 

k  =  R+  a*!®!  +  ai2^2  +  •••  +  \kSk 


where  K  <  M,  thereby  reducing  the  number  of  components  in  eigenvector  space  to 
estimate  Z^  by  M  -  K.  Values  of  K  and  of  2  and  3  can  produce  extremely 
good  reproductions  of  Z-,  while  greatly  reducing  the  number  of  parameters 
needed  to  describe  the  shape.  Further,  the  method  used  to  derive  the  parameters 
is  objective  and  rigorous. 

The  parameters  now  used  to  describe  the  shape  are  the  coefficients  aj,  a2, 
...,  a  ,  which  indicate  the  importance  of  the  shape  functions  e^,  e^,  ..., 
in  the  given  inlet  cross  section.  For  the  following  analysis,  M  varied  from 
20  to  60  and  N  varied  from  67  to  more  than  420.  In  the  cases  of  channel  cross 
sections  at  minimum  inlet  width  and  for  the  profiles  along  the  channel  thalweg, 
extra  profiles  were  analyzed  to  provide  a  wider  variety  of  inlet  conditions. 

It  should  be  noted  that  the  method  was  described  by  an  example  using  a  cross 
section  or  profile.  It  is  evident  that  for  any  single-valued  function  in  two- 
or  three-dimensional  space,  a  component  vector 

2=  [f(xlt  y  j) ,  f(x2,  y2),  f(xw,  yM)  ]  (8) 


can  be  defined  if  an  intrisic  grid  system  (x^,  y^)  can  be  established  for 
every  observation  of  f.  Hence,  as  is  the  case  for  the  ebb  tidal  delta  geom¬ 
etry  studied  here,  spatial  fields  of  data  can  be  analyzed. 

a.  Shape  of  the  Cross  Section  at  Minimum  Inlet  Width.  For  every  inlet 
examined,  the  depth  as  a  function  of  distance  from  the  left  side  of  the  minimum 
width  cross  section  was  recorded.  The  depth  d^  was  interpolated  linearly  at 
60  evenly  spaced  points  across  the  inlet,  resulting  in  a  60-component  vector 

W  =  Cd i ,  d2,  •  •  • »  d6o)  (9) 


Plots  of  the  cross  sections  are  provided  in  Appendix  A.  The  depths  were  normal¬ 
ized  by  the  minimum  width  W  of  the  inlet  to  give  a  dimensionless  geometry 


*  _  dj  £2  .  d60 

-  '  w-’  W  ’  •  •  •  »  W 


(10) 
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An  eigenvector  analysis  was  performed  and  the  first  three  shapes  plotted  in 
Figure  5  along  with  the  mean  shape.  The  variance  associated  with  e_1  is  86 

percent  of  the  total,  8  percent  with  e^,  and  3  percent  with  e_^.  Thus,  three 

vectors  explained  97  percent  of  the  variance  and  the  remaining  57  shapes  were 
ignored. 

The  parameters  (a3,  a2,  and  a3)  describing  the  shape  of  the  cross  section 
will  be  noted  as  EMI,  EM2,  and  EM3  and  are  the  weighting,  dimensionless 
coefficients  for  shapes  e  ,  e_^,  and  e_  ,  respectively.  As  a  guide  to  the 

interpretations  of  the  values  of  EMI,  EM2,  and  EM3,  the  following  gener¬ 
alizations  are  made: 

(a)  Positive  EMI-cross  section  shallower  than  the  mean. 

(b)  Negative  EMI-cross  section  deeper  than  the  mean. 

(c)  Positive  EM2— left  asymmetric. 

(d)  Negative  EM2-right  asymmetric. 

(e)  Positive  EM3-single  channel. 

(f)  Negative  EM3-center  shoal  with  two  side  channels. 

Examples  of  four  inlets  that  represent  these  siv  variations  are  provided  in 
Figure  7.  It  should  be  noted  that  in  Figure  7  the  profiles  are  plotted  in  real 
depth  not  in  normalized  depth.  Hence,  although  Port  Royal  Sound  and  Dewees  have 
depths  of  about  40  to  50  feet,  their  relative  depths  are  vastly  different  be¬ 
cause  Port  Royal  Sound  is  10  times  the  width  of  Dewees.  Thus,  EMI  is  0.0590 
for  Port  Royal  and  -0.0810  for  Dewees. 

The  eigenvector  representations  appear  to  be  realistic  expressions  of  cross- 
sectional  shape  variability.  The  first  eigenvector  scales  the  shape  according 
to  its  shallowness  or  deepness  with  regard  to  the  mean.  The  second  eigenvector 
displays  the  asymmetry  of  the  inlet  channels.  The  third  eigenvector  displays  a 
tendency  toward  a  single  or  a  multiple  channel  system.  The  eigenvector  analyses 
not  only  provide  a  succinct  mathematical  representation  of  cross-section  shape, 
but  the  shapes  derived  closely  resemble  major  shape  variations  generally  recog¬ 
nized  as  important. 

It  should  be  noted  that  the  sign  of  EM2  is  a  directional  quantity.  For 
the  analyses  later  in  this  report  the  magnitude  of  EM2  is  used  without  the 
sign.  This  is  done  because  the  direction  of  asymmetry  is  unimportant  unless 
correlation  is  made  to  another  directional  quantity. 

b.  Shape  of  the  Main  Channel  Depth  Profile.  The  method  for  determining 
the  centerline  of  the  main  channel  was  discussed  previously.  As  with  the  MIWC, 
depths  were  recorded  as  a  function  of  arc  length  along  the  channel  centerline 
from  the  minimum  width  cross  section  to  the  crest  of  the  ebb  tidal  delta.  At 
20  evenly  spaced  points  the  depth  was  linearly  interpolated  and  a  profile  vector 
constructed 


P  =  (dj,  d2,  ...,  d20)  (11) 

The  first  two  calculated  eigenvectors  are  plotted  with  the  mean  in  Figure  6. 
The  first  eigenvector  explained  87  percent  of  the  total  variance,  the  second  8 
percent,  and  the  third  3  percent.  Only  the  first  two  eigenvectors  are  used  to 
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Figure  7.  Examples  of  inlet  cross  sections  showing  the  major  variations 
in  EMI,  EM2,  and  EM3.  Note  that  the  eigenvectors  were 
computed  for  the  normalized  depth. 

parameterize  channel  shape,  and  the  other  18  are  ignored.  Four  examples  of 
natural  profiles  for  different  eigenvector  weightings  are  given  in  Figure  8. 

The  parameters  used  are  denoted  as  EC1  and  EC2  and  the  respective 
weighting  coefficients  for  e_  and  e^.  When  EC1  is  a  positive  number  the 
result  is  a  channel  which  is  in  general  shallower  than  the  mean  of  all  channel 
when  it  is  a  negative  value,  the  result  is  the  opposite.  When  EC2  is  a  posi 
tive  value,  the  inner  part  of  the  channel  (closer  to  the  MIWC)  is  relatively 
deep  and  the  outer  part  is  shallow.  For  EC2  negative,  the  inner  part  of  the 
channel  is  shallow  and  the  outer  part  deep. 

Again,  the  patterns  represented  in  the  eigenvectors  are  variations  typical 
of  inlet  channel  profiles.  The  first  eigenvector  is  very  much  related  to  the 
total  depth  of  the  channel.  The  second  eigenvector  indicates  the  presence  of 
a  bar  near  one  end  of  the  channel  and  scour  near  the  other  end. 

c.  Shape  of  the  Ebb  Delta.  A  primary  difficulty  in  parameterization  of 
inlet  morphology  is  representation  of  the  offset  of  the  inlet  and  the  shape  of 
the  ebb  tidal  delta.  Several  approaches  were  tried  unsuccessfully.  In  all 
cases  definitional  ambiguities  occurred,  particularly  in  the  definition  of 
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offset.  The  set  of  parameters  tended  to  be  unwieldy  to  analyze  and  very  diffi¬ 
cult  to  interpret.  An  eigenvector  analysis  of  the  spatial  pattern  of  depths 
over  the  ebb  delta  was  attempted  and  the  results  afforded  a  simpler  set  of 
parameters. 
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Figure  8. 


(1)  Construction  and  Alinement  of  Ebb  Tidal  Delta  Grid.  A  semicircu¬ 
lar  grid  (Fig.  9)  was  constructed  to  compare  the  morphology  of  the  67  ebb  tidal 
deltas  at  the  same  size,  i.e.,  to  normalize  the  delta  geometry.  A  radius  of  5 
inches  was  chosen  as  a  working  size  to  represent  the  prototype  distance,  r, 
as  shown  in  Figure  10.  The  radii  were  spaced  at  10°  intervals,  and  the  con¬ 
centric  semicircles  were  separated  by  one-half  inch. 


The  alinement  of  the  grid  was  determined  by  the  trend  of  the  local  shore¬ 
line  and  the  location  of  the  inlet  minimum  width  midpoint.  The  midpoint  of  the 
base  line  of  the  grid  was  positioned  on  the  midpoint  of  the  inlet  minimum  width 
line,  and  the  base  line  of  the  grid  was  set  parallel  to  the  general  trend  of 
the  local  shoreline.  The  radii  were  directed  seaward,  and  the  recording  points 
were  located  over  parts  of  the  ebb  tidal  delta.  The  left  side  of  the  grid 
viewed  from  the  base  line  facing  seaward  remains  left  for  all  coastlines;  i.e., 
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Figure  9.  Grid  mesh  for  describing  ebb  delta  geometry.  Sampling  points  are  numbered. 
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Figure  10.  Measurement  of  ebb  delta  radius,  r,  and  alinement  of  ebb  delta  grid. 


for  the  U.S.  Atlantic  coast,  left  will  be  directed  generally  northward  and  for 
the  U.S.  Pacific  coast,  left  will  be  trending  southward.  After  the  depths  were 
recorded  the  grid  was  "flipped,"  as  necessary,  so  that  the  shallower  side  of 
the  inlet  was  always  on  the  right  side  of  the  grid.  This  was  done  because  the 
objective  of  the  analysis  in  which  the  data  are  used  is  to  define  the  geometry 
of  the  delta,  which  is  by  definition  nondirectional .  If  a  grid  point  fell  out¬ 
side  of  the  delta,  the  depth  was  recorded  as  well.  Plots  of  the  ebb  deltas 
used  are  given  in  normalized  distance  format  in  Appendix  A. 

Because  of  the  differing  ebb  delta  sizes  it  was  necessary  to  normalize  the 
ebb  deltas  so  that  they  were  geometrically  similar;  i.e.,  the  ebb  delta  sizes 
on  the  map  were  reduced  or  expanded  to  have  approximately  the  same  area.  This 
was  achieved  by  reducing  (or  expanding)  the  scale  of  the  chart  until  all  of 
the  ebb  deltas,  defined  as  all  points  interior  of  the  depth  contour  equivalent 
to  the  controlling  depth  (DCC),  were  just  contained  within  the  semicircle  de¬ 
fined  previously.  This  was,  in  general,  straightforward  but  there  were  a  few 
instances  such  as  Moss  Bay,  California,  where  the  shoreline  is  not  linear  but 
makes  a  90°  bend  near  the  inlet.  In  such  cases  the  alinement  of  the  grid  was. 
adjusted  to  be  as  consistent  as  possible  to  the  other  charts. 

For  a  number  of  inlets,  the  basic  charts  did  not  contain  sufficient  data. 
Depths  from  charts  closest  in  time  to  the  basic  chart  used  were  substituted. 

The  inherent  assumption  in  this  substitution  is  that  the  rate  of  change  of  ebb 
delta  shape  is  small  with  time.  This  should  be  recognized  as  a  possible  source 
of  random  error  in  comparison  to  the  other  variates. 

(2)  Selection  of  Grid  Recording  Points.  Within  the  grid  constructed 
for  normalizing  the  ebb  tidal  deltas,  there  are  191  radii  semicircle  intersec¬ 
tions.  Twenty-two  were  selected  to  represent  the  topography  of  the  ebb  tidal 
delta.  The  location  of  those  points  is  shown  in  the  sample  grid  (Fig.  9). 

The  amount  and  location  of  the  recording  points  were  not  intended  to  give  a 
precise  description  of  the  delta  but  to  give  a  consistent  recording  of  data 
for  similar  locations  on  various  deltas  that  have  been  normalized  by  size. 

The  number  of  points  is  in  part  limited  by  the  number  of  deltas  available  for 
eigenvector  analysis  which  requires  that  the  number  of  grid  points  be  substan¬ 
tially  less  than  the  number  of  inlets  (22  versus  67). 

(3)  Eigenvectors  of  the  Ebb  Tidal  Delta  Geometry.  To  reiterate,  the 
geometry  analyzed  represents  a  reorientation  so  that  the  primary  mass  of  the 
delta  is  to  the  right  side  of  the  grid.  Further,  the  depths  over  the  delta 
have  been  normalized  by  the  average  depth  of  the  MIWC.  Thus,  the  ebb  delta 
vector  is  given  by 

£  =  (hj,  h2 ,  ...»  h22)  (12) 

where 


with  the  depth  at  the  i^  intersection  on  the  grid.  The  eigenvector 

analysis  was  performed  on  the  correlation  matrix  rather  than  the  covariance 
matrix  to  prevent  the  depths  around  the  perimeter  of  the  grid  from  dominating 
the  analysis. 


wi-th  N  the  number  of  samples,  a^,  aj  the  standard  deviations  and  d^-  and 
j  the  mean  depths  at  locations  i  and  j  respectively,  and  d^  and  d^- 
the  depths  at  locations  i  and  j  for  inlet  k.  The  mean  geometry  and  the 
first  three  eigenvectors  are  presented  in  Figure  11.  Examples  of  inlets  cor¬ 
responding  to  different  weighting  are  provided  in  Figure  12. 

The  mean  geometry  exhibits  a  central  channel  and  a  shoal  mass  on  the  right 
that  is  distinctly  shallower  than  the  left  as  is  expected  by  definition.  The 
first  eigenvector  explains  36  percent  of  the  total  variance  and  represents  a 
deep  delta  (depths  deeper  than  the  mean  delta)  when  the  weighting  coefficient, 
EDI,  is  positively  valued;  when  negatively  weighted  the  depths  are  shallow 
indicating  a  relatively  well-developed  delta.  The  second  eigenvector  explains 
15  percent  of  the  variation.  When  its  weighting  function,  ED2,  is  positive, 
the  right-hand  shoal  area  is  made  even  shallower,  the  central  channel  is 
deepened,  and  the  left-hand  shoal  is  shallower  (but  not  to  the  same  degree  as 
the  right-hand  shoal).  When  ED2  is  negative  the  shoals  are  deeper  and  the 
channel  shallower.  The  third  eigenvector  explains  10  percent  of  the  variance, 
is  somewhat  more  complex,  and  appears  to  represent  a  finer  representation  of 
the  channel  location.  Only  EDI  and  ED2  are  used  in  the  analyses  because 
the  third  eigenvector  appears  to  be  a  smaller  scale  variation  that  is  less 
likely  to  be  related  to  the  other  descriptors  used  here.  EDI  and  ED2  are 
dimensionless. 


To  further  discuss  the  geometry  represented  by  the  parameters  EDI  and 
ED2,  it  should  be  noted  that  the  values  analyzed  are  relative  depths  (depth 
divided  by  DMA).  Thus,  EDI  represents  variation  in  ebb  delta  thickness  rela¬ 
tive  to  the  average  channel  depth.  For  EDI  positive,  the  delta  is  relatively 
deep;  i.e.,  the  differences  between  channel  depth  and  shoal  depths  are  less 
than  when  EDI  is  negative.  For  EDI  <  0,  the  shoals  are  higher  relative  to 
the  channel  depth,  thus  to  a  large  degree  EDI  indicates  how  incised  into  the 
shoal  the  channel  is.  ED2  to  a  certain  degree  as  well  indicates  not  only  an 
increased  (or  decreased)  asymmetry  of  the  shoal  mass,  but  how  incised  the 
channel  is.  When  ED2  is  positive  the  shoals  are  more  asymmetric  and  the 
channels  deeper.  When  negative,  the  reverse  is  true. 

A  comparison  of  the  variance  explained  by  the  first  three  eigenvectors  in 
the  ebb  tidal  delta  analysis  to  that  for  the  channel  (first  two  eigenvectors) 
and  MIWC  (first  three  eigenvectors)  analyses  indicates  that  only  61  percent  as 
opposed  to  95-  to  97-percent  variance  is  explained.  This  is  due  to  two  major 
reasons.  The  ebb  delta  analysis  is  over  a  two-dimensional  grid  compared  to 
one-dimensional  traverse;  thus,  there  is  an  additional  degree  of  variation 
possible.  Secondly,  the  results  tend  to  indicate  that  the  channel  and  minimum 
inlet  width  cross-sectional  shapes  are  highly  organized  but  that  the  ebb  delta 
shape  is  less  so.  In  the  first  two  cases,  the  shape  factors  are  likely  to  be 
strongly  related  to  tidal  currents  which  are  perhaps  an  order  of  magnitude 
larger  than  the  wave-induced  forces.  Over  the  ebb  tidal  delta,  however,  the 
difference  is  less;  thus,  the  shape  must  strongly  represent  the  interplay  of 
two  driving  forces. 
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Figure  11.  Contour  maps  of  mean  ebb  delta  geometry  and  the  first  three  eigenvectors  of  ebb  delta  geometry 
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> .  Distribution  of  Inlet  Parameters. 

The  arithmetic  mean  and  standard  deviation  of  each  of  the  13  variables  are 
given  in  Table  2.  These  variables  include  both  the  direct  geometric  parameters 
(W,  L,  DMX,  DMA,  DDC,  and  AED)  and  the  eigenvector  parameters  for  the  cross- 
section  shape  (EMI,  EM2 ,  and  EM3) ,  channel  profile  (EC1  and  EC2 ) ,  and  ebb  delta 
(EDI  and  ED2).  A  listing  of  the  data  by  inlet  is  provided  in  Appendix  B.  Histo¬ 
grams  of  the  values  are  provided  in  Appendix  C. 


Table  2.  Mean  and  standard  deviation  of  parameters. 


Variable 

and  No. 

Mean 

Standard  deviation 

DMX 

1 

-28.9851 

17.4230 

DMA 

2 

-11.4851 

6.90464 

W 

3 

3332.04 

4036.59 

DCC 

4 

-8.51493 

6.65547 

L 

5 

7111.99 

5692.89 

EMI 

0 

0.0085268 

0.074429 

EM2 

7 

0.0032447 

0.032041 

EM3 

8 

-0.000441 

0.023791 

ECl 

9 

14.7692 

49. 3612 

EC2 

10 

5 . 75  760 

12.9366 

EDI 

11 

0.000002985 

2.89729 

ED2 

12 

0.000001492 

1.58685 

AED 

13 

2.52133 

4.40776 

III.  RELATIONSHIPS  AMONG  THE  GEOMETRIC  PARAMETERS 
1.  Procedures. 


The  number  of  possible  combinations  of  13  variables,  as  well  as  dimension¬ 
less  and  bifunctional  relationships  that  can  be  formed  from  these  variables,  is 
large.  Three  guidelines  were  used  to  reduce  the  task.  First,  all  combinations 
of  pairs  of  the  13  parameters  were  examined.  Then,  relationships  to  the  area 
of  the  MI WC  were  investigated.  Finally,  various  hut  not  all  dimensionless  com¬ 
binations  were  considered.  The  combinations  selected  were  those  that  appeared 
both  logical  and  fruitful  on  the  basis  of  relationships  seen  on  other  plots. 

The  procedure  used  to  evaluate  possible  relationships  was  to  first  plot  the 
variables  concerned  in  nont ransformed  coordinates.  If  a  functional  relation¬ 
ship  appeared,  the  coordinates  were  transformed  by  appropriate  combination  of 
logarithm  transforms  if  the  relationships  appeared  nonlinear.  In  the  appropri¬ 
ate  coordinate  system  a  linear  regression  was  performed  to  statistically  fit 
the  curve,  estimate  the  degree  of  fit  through  the  coefficient  of  determination 
(R?),  and  produce  95-percent  confidence  bands.  The  methods  used  are  common 
statistical  techniques  and  are  discussed  in  a  number  of  textbooks  (Kendall  and 
Stuart,  1901;  Krumbein  and  Graybill,  1905;  Dixon  and  Massey,  1909).  Based  on 
this  regression  analysis,  an  i  ratio  (the  ratio  between  the  mean  squares  due 
to  the  regression  and  the  mean  sum  of  the  squared  deviations  not  explained  by 
the  regression)  was  computed.  Because  there  are  07  inlets  in  the  study  this 
ratio  must  be  greater  than  a  7.08  value  for  the  F  distribution  with  vj  =  1 
and  vi?  =  05  degrees  of  freedom  for  the  regression,  or  curve  fit,  to  be  sig¬ 
nificant  at  a  5-percent  level. 
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The  analyses  are  discussed  first  for  direct  relationships  between  individ¬ 
ual  variables.  Next,  relationships  to  the  cross-sectional  area  of  the  MIWC 
are  discussed.  Finally,  relationships  to  a  dimensionless  parameter  W/L  are 
discussed. 

2 .  Direct  Relationships  Between  Various  Parameters, 
a.  Strong  Relationships. 

(1)  DMX  versus  DMA.  Figure  13  provides  a  linear  plot  of  DMX  versus 
DMA.  The  parameters  appear  strongly  related  as  is  confirmed  in  the  curve  fit 
analysis.  The  statistically  derived  relationship  is 

DMA  =  -1.42  +  0.347  DMX  (14) 

noting  both  DMA  and  DMX  are  negatively  defined  (i.e.,  a  depth  of  17  feet  is 
recorded  as  -17).  The  F  ratio  of  the  curve  fitting  regression  is  215  and  R2 
is  76.8  percent;  both  are  extremely  significant. 
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Figure  13.  DMX  versus  DMA. 

It  was  not  unexpected  to  find  a  relationship  between  the  average  depth 
of  the  MIWC  and  the  maximum  depth  therein.  However,  the  strength  of  the  rela¬ 
tionship  is  greater  than  initially  expected.  It  would  appear  that  there  is  a 
higher  degree  of  coherence  in  the  form  of  the  cross  section  than  generally 
a  s  s  ume  d . 


(2)  DMX  versus  DCC.  Figure  14  provides  a  log  plot  of  these  variates 
(signs  deleted).  The  relationship  found  is 

IH.C  =  (1.5662  DMX0-78  (15) 
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Figure  14.  DMX  versus  DCC. 

The  relationship  has  an  F  ratio  of  89  which  is  significant.  The  coefficient 
of  determination  is  57.9  percent. 

The  relationship  between  the  maximum  depth  in  the  MIWC  and  the  depth  at  the 
crest  of  the  outer  bar  again  is  highly  significant  and  implies  a  consistent 
adjustment  of  the  channel  depth  profile  to  the  crest  of  the  outer  bar.  It 
should  be  noted  however  that  considerable  variability  still  remains. 

C 3 J  DMX  versus  L.  Figure  15  provides  a  linear  plot  of  L  as  a  func¬ 
tion  of  DMX.  The  relationship  between  DMX  and  L  is  less  significant  than 
those  previously  discussed  but  is  still  highly  significant  with  an  F  ratio 
of  60  and  R2  of  48  percent.  The  curve  fit  is  (with  DMX  negatively  defined) 

L  =  539  -  226.7  DMX  (16) 

As  expected,  a  relatively  long  channel  is  associated  with  a  relatively  deep 
inlet  throat. 

(4)  DMX  versus  HCl.  Figure  16  is  a  linear  plot  of  these  parameters. 

If  I:C1  can  be  consistently  predicted,  the  channel  profile  can  be  predicted 
reasonably  well  because  F.Cl  represents  87  percent  of  the  shape  variation  in 
the  profile.  The  F  ratio  of  346  is  the  highest  achieved  in  this  study,  as 
is  the  coefficient  of  determination  value  of  84  percent.  The  linear  relation¬ 
ship  is  (with  DMX  negatively  defined) 

FC1  =  90.1  +  2.6  DMX  (17) 

where  DMX  is  again  negative  valued.  It  can  be  seen  that  a  relatively  deep 
channel  corresponds  to  a  relatively  deep  inlet,  remembering  that  a  channel 
deeper  than  the  mean  has  a  negative  TCI. 
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(5)  EC1  versus  DMA,  DCC.  Given  the  relationships  between  DMA,  DCC, 
and  DMX,  good  relationships  to  EC1  are  expected.  The  equations,  F  ratios, 
and  R2  values  are 

EC1  =  81.9  +  5.8  DMA  (18) 

with  F  =  131  and  R2  =  66.9  percent,  and 

EC1  =  66.7  +  6.1  DCC  (19) 

with  F  =  136  and  R2  =  67.6  percent.  DMA  and  DCC  are  negative  valued.  Thus, 
a  relatively  deep  channel  corresponds  to  both  a  relatively  deep  inlet  and  a 
relatively  deep  bar  channel.  The  corresponding  plots  are  Figures  17  and  18. 

(6)  EC1  versus  L.  Figure  19  provides  a  linear  plot  of  EC1  versus 
L.  The  relationship  is  strong  with  F  at  89  and  R2  at  58  percent.  The 
equation  is 

EC1  =  61.7  -  0.0067L  (20) 

In  this  case,  lengthening  the  channel  increases  the  channel  depth  and  the  con¬ 
trolling  depth  at  the  ebb  tidal  delta. 

(7)  All D  versus  L.  As  Figure  20  shows,  the  area  of  the  ebb  delta  (AED) 
is  strongly  related  to  channel  length  and  in  statistical  significance  is  second 
only  to  that  betwefen  DMX  and  EC1.  The  F  ratio  is  327  and  R2  is  83.4  per 
cent.  The  relationship  is  nonlinear  and  given  by 

AED  =  3.9245  *  10' 7L171  (21) 

AED  is  measured  in  square  miles,  L  in  feet.  Since  the  ebb  delta  is  bounded 
on  the  offshore  side  by  the  DCC  contour  and  the  channel  extends  to  this  same 
contour,  it  was  expected  that  increasing  L  would  also  increase  AED. 

The  strong  relationship  between  L  and  AED  is  unexpected.  A  small  par¬ 
tial  correlation  was  initially  supposed  but  not  to  the  degree  found.  The  lim¬ 
its  of  the  delta  are  defined  in  this  study  by  the  contour  equivalent  to  DCC, 
the  depth  of  the  crest  of  the  outer  bar  in  the  channel.  Given  the  relation¬ 
ships  described  previously  in  combination  with  this  relationship,  it  is  evi¬ 
dent  that  there  is  a  strong  covariance  among  many  of  the  major  components  of 
the  inlet  geometry.  It  is  interesting  to  note  that  IV  and  to  a  lesser  degree 
DMA  are  not  related  to  the  other  parameters. 

b.  Weak  Relationships.  The  relationships  discussed  above  are  all  highly 
significant.  There  are  a  number  of  other  relationships  involving  DMX,  DMA, 

DCC,  and  EDI,  ED2,  EMI  and  ED2,  EC1  and  ED2 ,  and  AED  and  EC  1 ,  which  appear 
potentially  promising.  In  all  cases,  however,  additional  data  are  needed  to 
further  define  the  curve  and  confirm  a  functional  relationship.  For  this 
reason  curve  fitting  was  not  performed. 


Figure  18.  DCC  versus  E Cl. 


From  Figures  21,  22,  and  23  it  can  be  seen  that  the  ebb  tidal  delta  func¬ 
tion  ED2  (which  describes  the  degree  of  delta  asymmetry)  appears  related  to 
DMX,  DMA,  and  DCC  in  the  following  manner 


Ca 

ED2  =  -2.  -  Ci,  B  >  0 

yb 


(22) 


where  y  is  either  DMX,  DMA,  or  DCC  (positively  defined)  and  Cq,  Cj,  and  B 
are  positive  constants  that  must  be  determined  in  a  regression  analysis.  Thus, 
for  decreasing  depths  in  the  channels,  shoal  asymmetry  appears  to  become  some¬ 
what  more  prevalent.  Channel  length  L  follows  a  relationship  with  ED2 
equivalent,  in  form,  to  those  followed  by  DMX,  DMA,  and  DCC  (see  Fig.  24). 

Figure  25  indicates  a  relationship  between  the  shape  of  the  MIWC  given  by 
EMI  and  the  relative  depth  of  the  outer  bar  given  by  EDI.  The  form  of  the 
relationship  (with  Fq,  Fi,  F2  and  A  positive  constants)  is 


EDI  =  F0  +  F 1 (EMI  -  F2)A,  A  >  0 


(23) 


This  indicates  that  for  relatively  shallow  inlets  (widths  greater  compared  to 
depths;  positive  EMI),  the  ebb  tidal  delta  is  relatively  flat  (depths  large 
compared  to  DMA;  positive  EDI).  A  less  veil-developed  ebb  tidal  delta  is 
related  to  a  shallow  channel,  with  the  channel  not  as  well  incised  into  the 
shoal.  With  increasingly  negative  EMI,  the  channel  is  more  incised  into  the 
shoal . 

A  final  relationship  is  shown  in  Figure  26  between  AED  and  EC1.  The 
form  of  the  relationship  with  Go,  Gj,  G2,  and  J  positive  constants 


Gd 

EC1  =  - - -  -  G2,  J  >  0  (24) 

(Gx  -  AED)17 


Given  the  correlation  between  L  and  AED  and  L  and  EC1,  this  is  not 
unexpected. 

3.  Relationships  to  the  Cross-Sectional  Area  of  the  MIWC. 

a.  Strong  Relationships.  The  geometric  parameter  previously  shown  to  be 
most  important  to  inlet  hydraulics  is  the  minimum  cross-sectional  area.  Its 
relationship  to  the  tidal  prism  of  the  inlet  has  been  defined  by  O'Brien  (1931) 
and  Jarrett  (1976).  As  mentioned  previously,  the  MIWC  is,  in  most  instances, 
located  close  to  the  minimum  area  cioss  section.  Thus,  the  areas  should  be 
approximately  equivalent.  Since  W  and  DMA  did  not  appear  as  dominant 
factors  in  the  previous  analyses,  it  was  decided  to  consider  their  bifunc¬ 
tional  relationship  to  the  area  of  the  MIWC,  Ac 

Ae  =  W  x  DMA  (25) 
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Figure  21.  DMX  versus  FD2, 
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Figure  22.  DMA  versus  FD2. 
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Figure  23.  DCC  versus  ED2 . 
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Figure  24.  I.  versus  LD2. 


Log  plots  of  kQ  versus  DMX,  DCC,  L,  and  AED  appeared  highly  significant.  The 
functional  relationships,  associated  F  ratios  and  R2  values  are 

(1)  DMX  =  0.5479A°-38,  F  =  133,  R2  =  67.6  percent. 

(2)  DCC  =  0.2367A8-34,  F  =  66,  R2  =  50.2  percent. 

(3)  L  =  23.92A8-55,  F  =  238,  R2  =  77.5  percent. 

(4)  AED  =  3.  1480  *  10-5A3-04,  F  =  260,  R2  =  80.0  percent. 

The  plots  are  provided  in  Figures  27  to  30.  These  relationships  are  unques¬ 
tionably  significant  statistically  and  provide  valuable  insight  to  the  excel¬ 
lent  covariance  relationships  noted  previously.  The  relations  have  significant 
design  implications  and  provide  an  indication  of  the  great  degree  of  the  co¬ 
adjustment  of  the  inlet  geometry.  A  broader  discussion  is  given  at  the  end  of 
this  section. 

b.  Weak  Relationships.  Figure  31  indicates  a  relationship  between  ED2 
and  \a.  It  is  approximated  by 


C0 

ED2  =  -£■  -  Ci,  B  >  0  (26) 

R<2 

with  C0,  Cj,  and  B  positive  constants  unrelated  to  any  other  constants  previ- 
osly  defined. 

4.  Relationships  to  W/L. 

Of  the  numerous  dimensionless  relationships  tried,  the  parameter  W/L  was 
most  successful.  The  primary  relationship  found  is  to  the  dimensionless  param¬ 
eter  DMA/DCC  shown  in  Figure  32.  The  relationship  is  described  by 

=  0.9289  (W/L)'0-42  (27) 

The  F  ratio  is  30.8,  R2  is  32.1  percent.  The  regression  is  significant 
above  5  percent.  There  is  appreciable  scatter,  but  the  low  R2  value  is 
likely  due  to  three  to  five  points.  Additional  points  would  possibly  increase 
the  R2  value. 

There  is  a  weak  relationship  (Fig.  33)  of  W/L  to  EDI  (which  is  dimen¬ 
sionless).  It  is  given  by 

EDI  =  -1.  36  +  2. 85  (W/L)  (28) 

The  F  ratio  is  8.4  and  just  barely  significant. 

Both  relationships  tend  to  indicate  the  following  adjustment.  As  W/L 
decreases  from  a  value  of  1.5  or  so,  the  ebb  tidal  delta  tends  to  become  more 
developed  (EDI,  smaller)  and  a  strong  bar  crest  develops  (DMA/DCC  increases). 
This  is  shown  graphically  in  Figure  34,  which  provides  average  channel  profiles 
with  the  depths  normalized  by  DMA,  grouped  -by  classes  of  W/L  values. 
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5.  1)  is  cuss  ion. 


The  analyses  presented  show  strong  relationships  among  many  of  the  parame¬ 
ters  selected  to  describe  inlet  geometry.  Weak-to-moderate  partial  correlations 
had  been  expected,  but  the  strength  of  the  relationships  implies  a  more  highly 
organized  covariance  or  coadjustment  of  many  elements  of  inlet  geometry  than  has 
been  recognized  previously.  The  implications  of  the  results  have  inpact  upon 
the  understanding  of  the  adjustment  of  inlets  to  wave  and  tide  processes  and 
may  indicate  relationships  useful  for  the  design  of  inlet  improvements . 

The  relationships  among  DMX,  DMA,  DCC,  L,  and  TCI  and  ATD  are  placed  into 
prospective  by  the  relationship  of  these  variables  to  the  area  of  the  MIWC, 

A,.  To  reiterate,  the  cross-sectional  area  is  directly  related  to  the  channel 
length,  depth  at  the  crest  of  the  outer  bar  in  the  channel,  and  t lie  ebb  tidal 
delta  area.  The  relationship  of  all  these  variables  suggests  that,  as  in  the 
case  of  the  tidal  prism,  the  cross-sectional  area  of  the  inlet  appears  to  be  a 
controlling  variable.  Thus,  for  a  given  tidal  prism,  the  area  is  to  a  certain 
degree  determined  for  fixed  tidal  range.  With  increasing  discharge,  kQ  in¬ 
creases,  and  as  long  as  velocities  remain  above  the  critical  velocity  for  sedi¬ 
ment  transport,  the  flow  in  the  channel  tends  on  the  average  to  be  large  enough' 
to  control  the  length  of  the  channel,  and,  therefore,  the  location  of  the  crest 
of  the  outer  bar.  The  strength  of  the  relationship  between  DCC  and  Aa  would 
indicate  that  even  at  this  point  where  wave  and  tidal  transport  would  concep¬ 
tually  be  of  the  same  order,  the  tidal  forces  (represented  by  Aa)  are  still 
predominant  in  influencing  certain  geometric  parameters.  Using  the  87-percent 
explanation  of  channel  profile  variance  by  TCI,  t he  control  of  the  channel 
deptli  profile  (represented  by  TCI,  TC2)  is  primarily  related  to  Ac  (given 
the  strong  relationship  between  A(.  and  TCI).  TC2  is  not  related  to  Aa  and 
it  is  suspected,  given  the  shape  of  TC2  (Tig.  b),  this  variable  represents 
an  adjustment  of  the  channel  for  varying  wave  climates.  However,  since  there 
is  no  parameter  distinctly  representing  wave  conditions  in  this  study  (which 
would  correspond  to  A„  for  tides),  it  is  impossible  to  further  evaluate  a  wind- 
wave  dependence. 

The  strength  of  relationship  between  A.  and  ATI)  is  surprising  even 
though  partial  correlation  is  expected.  However,  given  the  adjustment  of 
channel  length  to  A,,  and  the  tendency  for  the  principal  ebb  flow  to  remain 
as  a  jet,  simple  constraints  of  geometry,  wave  refraction,  and  diffraction 
appeared  to  dictate  the  eventual  size  of  the  bar,  as  discussed  by  Bates  (1953). 
Deviations  from  this  gross-scale  geometry  would  appear  to  result  from  major 
geologic  differences  such  as  deviation  away  from  a  fairly  straight  coastline. 

The  correlation  between  the  second  eigenvector  of  ebb  tidal  delta  geometry 
(TDJj  and  the  cross-sectional  area  shows  that  for  smaller  cross-sectional 
areas  the  right  side  (in  normalized  geometry  space)  of  ebb  tidal  delta  is  rela¬ 
tively  more  highly  developed  than  for  larger  A,,.  This  would  agree  with  the 
scaling  of  the  tide  and  wave  processes  discussed  before.  For  the  smaller  , 
the  discharges  are  lower,  and  velocities  over  the  inlet  delta  flanks  are  proba¬ 
bly  less.  Thus,  for  even  moderate  wave  conditions,  sediment  transported  onto 
the  shoals  is  not  moved  as  easily  offshore  or  to  edge  of  the  delta. 

The  analysis  shows  that  the  maximum  depth  in  the  cross  section  is  highly 
related  to  A,.  This  would  seem  to  imply  that  in  the  principal  flow  area  of 
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the  channel  the  geometry  of  the  channel  is  largely  determined  by  the  magnitude 
of  the  flow,  since  ka  has  been  shown  to  be  related  to  tidal  prism.  Although 
considerable  variability  in  the  shape  is  expected,  the  gross-scale  geometry  of 
the  channel  represented  mainly  by  DMX  appears  determined. 

It  is  interesting  to  note  again  that  W  and  (to  a  lesser  degree)  DMA  did 
not  exhibit  strong  covariance  or  relate  as  well  to  the  other  parameters  as  some 
of  the  other  parameters  did.  This  suggests  that  W  and  DMA  are  free  param¬ 
eters  with  respect  to  the  tide;  i.e.,  given  a  particular  cross  section,  W  and 
DMA  are  much  freer  to  adjust  themselves  and  may  widely  vary  according  to  wave 
climates.  It  is  unclear,  however,  how  DMA,  DMX,  kQ,  or  W  coadjusts  to  the  wave- 
tide  regime  as  it  is  unclear  how  and  the  tidal  prism  coadjust. 

The  preceding  comments  must  also  consider  that  the  cross-sectional  area 
A^  is  an  adjustment  to  capacity  of  the  tidal  current  and  wave  transport  of 
sediment  although  it  is  likely  that  for  an  kc  of  any  appreciable  size  the 
tidal  capacity  must  dominate.  The  results  presented  here  suggest  that  to  the 
first  order  the  tidal  flows,  scaled  by  ka,  determine  the  gross-scale  geometry 
of  the  inlet  delta  and  channels  (within  geologic  limitations).  The  wave-related 
changes  appear  as  a  modification  to  this  geometry  until  ka  becomes  so  small 
that  ebb  flows  are  weak  relative  to  the  wave  forces. 

The  relationships  developed  here  appear  to  have  potential  use  in  engineer¬ 
ing  planning  and  design.  They  represent  geometric  adjustments  to  natural  con¬ 
ditions,  and  at  present,  care  must  be  taken  in  their  use  until  the  implications 
suggested  by  the  relationships  have  been  thoroughly  analyzed. 

IV.  CLASSIFICATION  OF  INLETS 
1.  Mathematical  Considerations. 

Classification  is  essentially  a  statistical  process  whether  done  numeri¬ 
cally  or  manually.  If  the  set  of  all  possible  inlet  morphologies  (of  which 
this  study  has  a  small  sample  of  individuals)  is  considered,  the  hypothesis 
arises  that  there  are  K  subpopulations  of  inlets  which  have  different  mor¬ 
phologic  characteristics  due  to  some  basic  physical  difference  in  the  inlets 
themselves.  Two  problems  need  consideration: 

(a)  Can  the  characteristics  of  the  subpopulations  be  estimated 

through  proper  parameter  selection  and  sampling? 

(b)  Can  each  individual  in  the  sample  be  assigned  to  the  correct 

subpopulation? 

The  size  of  K  is  not  known,  nor  is  a  sufficient  set  of  parameters,  Xj,  X2, 

...,  Xn,  known  which  is  necessary  for  solution  of  the  problem.  Compounding 
these  basic  difficulties  are  random,  hopefully  unbiased,  errors  in  estimates 
of  Xj ,  X2 ,  ....  Xn  for  each  individual  in  this  sample  which  may  cause  mis- 
assignment  of  the  individual.  Finally,  it  must  be  assumed  that  the  sample 
used  contains  enough  examples  from  the  K  subpopulations  to  allow  resolution 
of  the  problem. 

Because  of  these  difficulties,  objective  statistical  techniques  for  classi¬ 
fication  analyses  were  only  recently  developed,  with  the  biologic  scientific 
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community  in  the  forefront  of  its  development.  The  analysis  of  the  structure 
of  populations  through  objective  methods  is  termed  numerical  taxonomy.  Classi¬ 
fication  analyses  denote  derivation  of  taxonoroetric  structure.  Discrimination 
analyses  denote  assignment  of  individuals  to  the  classes  derived.  Sneath  and 
Sokal  (1973)  provide  an  excellent  introduction  to  the  subject;  a  simple  intro¬ 
duction  is  given  in  Davis  (1973). 

Numerical  taxonomists  have  developed  their  own  terminology  which  in  general 
will  not  be  applied  here.  The  objective  of  this  part  of  the  report  is  to  pre¬ 
sent  as  simply  as  possible  the  application  of  cluster  analysis  to  the  inlet 
classification  project.  The  reason  for  use  of  this  technique  to  help  unravel 
the  taxonometric  structure  of  inlet  morphologies  lies  in  he  sheer  mass  of 
information  that  must  be  analyzed  and  in  a  desire  to  be  ODjective  in  the  final 
stratification  of  inlets. 

2 .  Cluster  Analysis. 

The  method  employed  here  is  a  weighted  pair  group  average  (WPGA)  cluster¬ 
ing  technique  (Sneath  and  Sokal,  1973;  Davis,  1973).  If  N  parameters  are 
measured  for  M  inlets  the  data  can  be  represented  as  an  N  x  M  matrix  F 
in  which  the  element  is  the  measurement  of  the  j th  parameter  for  the 

i^n  inlet.  The  mean  value  of  each  parameter  can  be  calculated 

-  1  M 

fJ  =  M  ,  1  fkJ  C29) 

k=  l 

and  the  associated  standard  deviation  estimated 


=  M 


A  matrix  F°  (with  elements  f\-  •)  can  be  formed  by  transforming  each  element  of 
F  in  the  following  way  d 


f  •  ■  -  f  • 

fC  .  = 


F°  represents  transformation  of  parameters  to  have  mean  0  and  standard  devia¬ 
tion  1.  Finally,  a  matrix  C  can  be  computed  with  elements  C £j  defined  as 


w  ‘  s(j,  «?<=  - 


C^j  represents  an  buclidean  distance  between  inlets  i  and  j  in  the  M- 
dimensional  space  defined  by  the  M  parameters  standardized  to  have  mean  0 
and  standard  deviation  1.  Other  distance  measures  such  as  correlation  coeffi¬ 
cients  can  be  used.  Here  the  distance  function  is  chosen  because  its  interpre¬ 
tation  is  simple  and  geometrically  appealing.  In  most  cases  taxonometric 
structures  derived  using  the  distance  and  correlation  functions  are  equivalent 
(Sneath  and  Sokal,  1973;  Davis,  1973). 


Cluster  analysis  is  so  termed  because  the  technique  orders  the  individual 
inlets  into  groups  or  clusters  based  upon  certain  fixed  methodologies.  As 
indicated  in  Sneath  and  Sokal  (1973)  there  is  a  wealth  of  possibilities  from 
which  to  choose.  In  the  WPGA  method,  the  matrix  C  is  surveyed  to  find  the 
element  C^j  which  is  minimum.  The  other  elements  of  C  involving  either 
i  or  j  (such  as  or  C j0 ,  are  replaced  with  1/2  +  Cj^).  In  essence 

the  individual  inlets  i  and  j  are  replaced  by  a  composite,  synthetic  inlet 
equal  to  their  average.  The  inlets  joined  together  are  those  "closest"  or  most 
similar  in  terms  of  the  distance  function.  Tlius,  the  inlets  i  and  j  form 
a  cluster.  This  process  is  continued  with  new  clusters  being  formed  or  with 
inlets  added  to  old  clusters.  eventual  1 . ,  these  groups  of  inlets  are  likewise 
joined  together  until  the  distance  relationship  among  all  inlets,  as  arranged 
in  hierarchial  order  of  similarity,  has  been  determined.  The  pattern  in  which 
the  inlets  cluster,  and  the  values  of  the  function  C/J  at  which  the  cluster¬ 
ing  occurs,  can  be  shown  by  a  dendrogram  (Fig.  35). 

The  value  at  which  the  clustering  occurs  is  indicative  of  the  simi¬ 

larity  of  the  two  elements  under  consideration.  For  Cfj.  near  zero,  the  in¬ 
lets  are  quite  similar  because  the  distance  in  the  normalized  parameter  space 
is  small.  The  unresolved  problem  of  many  cluster  analyses  is  to  determine  at 
what  values  of  distinct  cluster  discrimination  occurs.  In  general,  the 

question  cannot  be  answered.  In  the  case  of  M  normally  distributed  random 
variables  with  mean  0  and  standard  deviation  1,  a  mean  value  of  the  Euclid¬ 
ean  distance  measure  C ij  can  be  calculated  for  which  the  hypothesis  that  in¬ 
dividual  clusters  are  significantly  different  is  acceptable.  In  the  case  of 
this  study  (for  M  =  13  variables),  this  value  is  approximately  1.3  with  a  95- 
percent  confidence  band  of  0.95  to  1.85  (Sneath  and  Sokal,  1973).  The  hypoth¬ 
esis  of  normal  distribution  certainly  is  not  true  for  all  of  the  variables 
treated  here,  but  the  precise  effect  of  this  upon  the  distance  function  is 
unknown.  Because  the  variables  do  not  differ  radically  from  the  normal,  the 
effect  is  not  expected  to  be  major.  However,  the  cluster  analysis  should  be 
considered  only  as  a  guide  to  possible  subpopulations. 

The  dendrogram  (or  graphical  presentations  of  the  taxonometric  structure) 
of  the  inlet  relationships  was  examined.  Dendrograms  are  constructed  so  that 
it  is  possible  to  determine  how  the  individual  inlets  are  grouped  together  as 
clusters,  at  what  values  of  the  distance  function  the  inlets  join  a  cluster 
and  the  values  of  the  distance  function  between  clusters.  Clusters ’with  dis¬ 
tance  function  values  of  1.3  or  greater  were  separated.  For  each  pair  of 
clusters,  the  hypothesis  that  the  difference  in  means  is  zero  was  tested  for 
each  of  the  13  variables  using  a  standard  student's  t  test  applied  for  compar¬ 
ison  of  means  (Davis,  1973).  Rejection  of  this  hypothesis  at  a  10-percent  level 
of  significance  will  be  termed  significant  clustering  on  the  basis  of  the  vari¬ 
able  tested. 

3.  Inlet  Clusters  Derived. 

The  dendrogram  showing  the  taxonometric  structure  of  the  67  inlets  (Fig.  35) 
is  explained  as  follows:  The  horizontal  axis  labeled  distance  coefficient  indi¬ 
cates  the  value  of  C^j  for  which  an  inlet  joins  one  cluster  or  two  clusters 
join  together.  The  vertical  axis  contains  the  names  of  the  inlets  and  the  values 
at  which  the  inlet  joins  the  dendrogram.  The  dotted  lines  are  the  branches  of 
the  dendrogram  displaying  how  the  inlets  are  linked.  In  Figure  35,  Moriches 
joins  Stump  Inlet  at  a  value  of  0.27,  and  Price  joins  this  cluster  with  a  value 
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35.  Dendrogram  of  the  cluster  analysis. 
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of  0.41.  The  cluster  containing  Moriches,  Stump,  Price,  Gordan,  and  Midnight 
is  joined  with  Tubbs  and  Fort  George.  The  distance  values  are  all  below  1.3, 
which  is  the  criterion  selected  as  the  point  of  inlet  cluster  definition.  The 
cluster  labeled  cluster  1  joins  cluster  2  with  a  value  of  1.34.  Because  the 
distance  function  is  above  1.3,  tiie  two  clusters  are  considered  far  enough 
apart  to  lie  statistically  different,  thus  indicating  two  subpopulations.  For 
the  distance  value  of  1.3  chosen  to  indicate  significant  clustering,  six  clus¬ 
ters  were  defined.  The  inlets  in  each  cluster  are  listed  in  Table  3.  Of  67 
inlets  in  the  analysis,  all  but  five  are  in  well-defined  clusters.  These  five 
inlets  will  be  termed  outliers  because  they  do  not  fit  in  the  clusters  defined. 
The  five  will  be  discussed  at  the  end  of  the  section  describing  the 

individual  clusters. 


Table  3.  Inlets  by  cluster  group. 


Croup  1 

Group  2 

Group  3 

Group  4 

Group  S 

Group  6 

Outl iers 

Moriches 

Stump 

Price 

Go  r Jon 

Midnight 

Tubbs 

Port  George 
Shal lotte 
Schast ian 
Netarts 

Corson 

Gaspari 11a 
North 

South  Santee 
Ponce  do  Leon 
Si  us  1  aw 
Clearwater 
Chincotoague 
Pass  A  Grille 
San  Luis 
ilattcras 

Big  Sarasota 

P i re  I s 1  and 

P>ri  gantine 
Oregon 

Hereford 

Great  Egg 

Coos  Bay 

Beach  Haven 
Stono 

Nassau 

Pnpps 

St.  Augustine 
Captiva 

Umpqua 

Yaquina 

Murrells 

Carolina  Beach 
Bodega 

Boca  Raton 
Bolinas 

Lockwoods  Folly 
Morro  Bay 

Alsea 

Nehalem 

! 

Townsend 

Metomkin 

Deuces 

Lighthouse 

Big  Marco 
Longboat 
Little  River 
Capers 

Drakes 

Redfish 

Wachapreague 

Siletz 

Beaufort 

Boca  Grande 
Pensacola 
Humboldt 
Tillamook 

Hi  1 lsboro 
Doboy 

Pass  Cavallo 
Willapa 

Grays  Harbor 

a.  Statistical  Measures  of  the  Significance  of  the  Cluster  Analysis.  The 
clusters  were  selected  on  the  basis  of  a  Euclidean  distance  value  of  1.3  which 
was  chosen  on  the  basis  of  the  expected  difference  if  the  variables  conformed 
to  the  constraints  listed  previously. 

Given  six  clusters  there  are  15  possible  pairs  of  clusters,  i.e.,  clusters 
1  and  2,  1  and  3,  and  so  forth.  In  Table  4,  for  each  of  the  15  pairs,  the  lev¬ 
els  of  significance  based  on  a  student's  t  test  for  the  hypothises  that  the  mean 
cluster  value  for  each  of  the  13  parameters  is  different  are  provided.  The  num¬ 
ber  of  times  that  this  hypothesis  is  accepted  as  true  provides  a  measure  of  how 
different  the  two  clusters  in  a  pair  are.  If  only  a  few  of  the  13  parameters 
test  out  as  significant,  then  statistically  there  is  less  confidence  that  the 
clusters  are  distinct.  Of  the  15  pairs  of  clusters,  clusters  1  and  3  are  most 
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Table  4.  Significance  of  differences  among  clusters,  based  on  a 


student’s  t  test  for  differences  between  means. 

riiKt.r  n:i  i  rs  ^ 


Variable 

and  No. 

Cluster 

pairs^ 

| 

1-2 

1-3 

1-4 

1-5 

ESI 

ESI 
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ca 
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■a 
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1 

5 
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X 

1 

1 

1 

X 

X 

1 

5 

i 

BS 

10 

1 

1 

DMA 

2 

1 

wm 

1 

1 

1 

1 

X 

X 

X 

1 

i 

X 

5 

5 

w 

3 

X 

mm 

S 

1 

X 

1 

1 

1 

10 

X 

5 

S 

S 

1 

DCC 

4 

X 

i 

X 

X 

1 

5 

X 

X 

1 

5 

1 

10 

1 

1 

L 

S 

5 

i 

5 

X 

1 

5 

S 

X 

X 

X 

1 

s 

1 

1 

I. Ml 

6 

1 

i 

1 

1 

S 

X 

1 

1 

X 

10 

5 

X 

1 

1 

/I.M2/ 

7 

10 

i 

1 

1 

; 

5 

1 

X 

10 

1 

X 

1 

1 

X 
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8 

i 

10 

i 

1 

X 

10 

1 

5 

S 

1 

10 

1 

5 

1 

fcCl 

9 

10 

s 

X 

10 

1 

BS 

X 

1 

5 

1 

1 

X 

1 

1 
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10 

X 

X 

1 

D 

BS 

BS 

S 

X 

10 

5 

X 

X 

X 

LD1 

11 

10 

1 

1 

BS 

10 

10 

10 

X 

mm 

5 

X 

tD  2 

12 

1 

X 

X 

KS 

BS 

5 

X 

5 

1 

BS 

1 

X 

AfcD 

13 

BS 

s 

X 

X 

mm 

5 
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10 

10 

X 

5 

BS 

10 

1 
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BS 

Mm 

mm 

SB 

10 

10 

6 

7 

BS 

19 

19 

BS 

n 

12 

9 

At  10 
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BS 

19 

ss 

ss 

l 

1 

0 

2 

BS 

19 

BS 

Bfl 

BB 

BB 

KM 

At  S 

pet 

Bfl 

B8 

ss 

BS 

MS 

S 

2 

2 

Bfl 

mm 

BS 

BS 

Bfl 

Bfl 

Bfl 

At  l 

pet 

MS 

SH 

bs 

;SS 

BS 

4 

4 

3 

.BS 

BS 

.BS 

BB 

.SB 

BS 

.BB1 

*  ■  no  difference 
10  ■  significant  above  10  pet 
S  ■  significant  above  5  pet 
1  •  significant  above  1  pet 


different  in  that  12  of  13  variables  have  mean  values  accepted  as  different  at 
a  10-percent  level  of  significance;  only  5  variables  are  significant  above  1 
percent,  however.  Clusters  1  and  5  have  eight  variables  different  above  1 
percent.  Comparisons  of  clusters  1  and  6,  3  and  5,  4  and  6,  and  S  and  6  have 
seven  variables  significantly  different  above  1  percent.  Clusters  2  and  4 
appear  most  similar  with  only  six  variables  different  at  significant  levels 
above  10  percent;  however,  four  of  the  variables  are  significantly  above  1 
percent. 

The  comparisons  presented  in  Table  4  strengthen  the  conclusion  that  the 
clusters  selected  stratify  the  inlet  sample  into  apparent  distinct  groups  which 
have  measurable  differences.  It  is  an  inability  of  current  statistical  theory, 
however,  to  state  with  a  level  of  certainty  that  the  groups  chosey  are  truly 
optimal  in  a  population-wide  sense.  However,  the  grouping  appears  a  useful 
separation  of  inlets  into  groups  which  are  distinguishable  in  terms  of 
geometry. 

An  interesting  observation  is  apparent  in  the  analysis  shown  in  Table  4. 
There  is  no  single  variate  that  serves  as  a  distinguishing  parameter  for  all 
clusters.  The  differences  between  clusters  result  from  differing  combina¬ 
tions  of  the  13  variates.  This  observation  agrees  in  large  degree  with  the 
observations  of  Sneath  and  Sokal  (1973),  who  note  that  it  is  a  rarity  to  find 
one  single  parameter  which  serves  as  the  distinguishing  character  in  popula¬ 
tions  characterized  by  continuous  variates.  This  condition  appears  to  be  a 
property  of  natural  systems  which  are  in  some  sense  evolutionary. 

The  primary  reasons  for  performing  the  cluster  analysis  were  to  achieve  an 
objective  classification  analysis  and  to  have  an  automated  analysis  because  of 
the  difficulty  of  trying  to  analyze  manually  many  variables  for  a  large  sample. 
The  problem  of  the  massiveness  of  the  data  set  and  its  multivariate  nature, 
though  circumvented  in  the  cluster  analysis,  is  still  present  when  the  charac¬ 
teristics  of  the  clusters  must  be  presented.  Table  5  provides  a  summary  of  the 
mean  value  and  standard  deviation  of  each  variable  by  cluster  group.  Figure  36 
provides  plots  of  these  values  by  group.  Figure  37  provides  plots  of  the  limits 
of  each  cluster  in  selected  bivariate  spaces.  Table  5  and  Figures  36  and  37  are 
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Figure  36.  Mean  and  standard  deviations  of  the 
13  variables  by  clusters. 


intended  to  help  clarify  differences  and  similarities  between  clusters.  The 
limits  are  defined  to  be  those  rectangles  with  centers  located  at 

where  is  the  mean  value  of  the  b ^  (of  13)  variates  for  the  a^ 

(of  6)  clusters.  The  comers  of  the  rectangles  are  located  at 

±  *  k^i  - 

where  aS^,  is  the  standard  deviation  associated  with  aX^.  Thus,  for  Figure 
37(a)  the  six  rectangles  corresponding  to  the  six  clusters  are  plotted  in  (W,L) 
space,  bach  rectangle  can  be  interpreted  as  the  region  in  (W,L)  space  for 
which  no  less  than  66  percent  of  the  inlets  in  each  cluster  lie.  Where  two 
rectangles  such  as  for  clusters  1  and  3  do  not  overlap,  it  is  evident  that 
based  on  W  and  L  the  clusters  are  relatively  disjointed,  for  the  same 
value  range  for  a  variable,  as  shown  in  Table  4.  Conversely,  overlapping 
rectangles  usually  indicate  that  there  is  no  significant  difference  between 
the  two  variables  (e.g.,  see  N  and  L  for  cluster  pairs  3  and  4  in  Table  4). 

The  degree  of  overlap  determines  the  degree  of  similarity  between  the  two 
clusters.  It  should  be  noted  that  major  overlap  can  occur  only  in  one  vari¬ 
able,  as  in  clusters  1  and  6  for  channel  length,  L.  Further,  even  though 
clusters  2  and  5  do  not  overlap  for  W,  they  do  in  L. 

b.  Description  of  the  Inlet  Clusters  (mean  values  are  given  for  all 
variates) . 


(1)  Cluster  1.  The  first  cluster  is  by  far  the  largest  and  represents 
almost  half  of  the  inlets.  The  geographic  extent  of  the  31  inlets  includes  in¬ 
lets  on  all  three  coasts.  The  dendrogram  (Fig.  35)  does  not  indicate  a  neces¬ 
sarily  homogeneous  grouping  because  of  the  number  of  small  clusters.  However, 
the  low  values  of  the  distance  coefficients  suggest  that  this  internal  strati¬ 
fication  is  possibly  due  to  only  a  few  variates.  Since  it  is  a  substructure, 
it  will  not  be  further  considered  here  because  it  occurred  at  too  low  a  level 
of  signi ficancc. 

From  Table  5,  a  cluster  1  inlet  has  moderate  width  and  length  (W  -  3,800 
feet,  L  -  5,800  feet).  In  (DMX, DMA)  space  (Fig.  37, b)  it  has  shallow  average 
and  maximum  depths  (DMA  -  7.5  feet,  DMX  -  -22  feet).  Depth  at  the  crest  of  the 
bar  is  moderate  (DCC  -  -7  feet).  The  MIWC  shape  is  quite  shallow  (FM1  -  0.053), 
tends  not  to  be  overly  asymmetric  ( 12112  =  0.0074),  and  is  neither  consistently 
deep-centered  nor  mul  t  ichanncled  (liM3  -  0.0018).  The  channel  profile  (Fig. 

37, f)  is  moderately  shallow  (LCl  -  31)  and  has  a  slightly  steeper  slope  (EC2 
-  3.4).  The  ebb  delta  geometry  is  moderately  deeper  (FD1  -  0.60)  and  not  ex¬ 
tremely  asymmetric  (ED2  ~  0.075).  The  ebb  delta  area  is  small  (Alii)  ~  1.45 
square  miles).  The  variables  by  which  cluster  1  differs  from  the  other  clusters 
are  given  in  Table  4. 

(2)  Cluster  2.  The  second  cluster  is  comprised  of  six  inlets  located 
on  the  three  coasts.  Cluster  2  inlets  are  moderately  wide  (W  =  2,400  feet), 
but  not  as  wide  as  cluster  1  inlets,  and  have  a  fairly  long  channel  (L  =  8,500 
feet).  The  MIWC  is  moderately  deep  (DMA  =  16.5  feet,  DMX  =  32  feet).  The 
depth  of  the  crest  of  the  ebb  tidal  delta  is  moderate  (DCC  =  8  feet).  File 
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MIWC  is  near  the  mean  shallowness  (LM1  =  D.009),  and  is  moderately  asymmetric 
(EM2  =  0.013,  1:M3  =  -0.031J.  The  channel  profile  is  relatively  shallow  but 
close  to  the  mean  (EL'l  =  13. 3J.  The  channel  profile  steepness  is  relatively 
variable  (LC2  =  0.59).  The  ebb  tidal  delta  is  moderately  deep  (bill  =  0.93J 
and  not  highly  asymmetric  (LD2  =  -1.5).  The  average  ebb  tidal  delta  area  is 
intermediate  (AED  =  1.98  square  milesj. 

(3)  Cluster  3.  Cluster  3  has  only  four  inlets.  Two  are  on  the 
Atlantic  coast  and  two  on  the  Pacific  coast.  These  inlets  have  the  narrowest 
width  (W  =  480  feet)  and  the  shortest  channel  length  (L  =  1,700  feet).  The 
depth  values  are  shallowest  (DMX  =  11  feet,  DMA  =  5  feet,  DCC  =  3  feet);  the 
channel  profile  is  also  shallow  (EC1  =  64,  EC2  =  -2.4).  The  MIWC  shape  is 
slightly  deeper  than  the  mean  shape  (EMI  =  -0.013)  and  more  asymmetric  (EM2  = 
0.036,  F.M3  =  -0.008).  The  ebb  tidal  delta  is  shallow  (EDI  =  -0.94),  fairly 
asymmetric  (ED2  =  3.35),  and  small  (AED  =  0.09  square  miles). 

(4)  Cluster  4.  Cluster  4  is  comprised  of  four  inlets,  three  of  which 
are  on  the  Pacific  coast.  The  minimum  inlet  width  is  relatively  small  (W  =  625 
feet)  and  the  channel  length  is  short  (L  =  2,850  feet).  The  principal  depths 
are  generally  in  the  slial low-to-intermediate  range  (DMX  =  25  feet,  DMA  =  11 
feet,  DCC  =  8  feet)  and  the  channel  profile  is  in  the  intermediate  value  ranges 
(EC1  =  26.4,  EC2  =  8.8).  The  MIWC  shape  is  relatively  deep  (EMI  =  -0.10)  and 
very  asymmetric  (HM2  =  0.099,  EM3  =  -0.045).  The  ebb  tidal  delta  is  relatively 
shallow  and  slightly  asymmetric  (EDI  =  -2.95,  ED2  =  0.17),  and  has  an  area  of 
1.18  square  miles. 

(5)  Cluster  5.  Cluster  5  includes  12  inlets  that  are  geographically 
distributed  on  all  three  coasts.  These  inlets  are  relatively  narrow  but  have 
moderately  long  channels  (W  =  990  feet,  L  =  6,300  feet).  Depths  in  the  channels 
are  relatively  deep  (DMX  =  32  feet,  DMA  =  15  feet,  DCC  =  6  feet,  EC1  =  18.5, 

EC2  =  13.6).  The  MIWC  geometry  is  deep  and  moderately  asymmetric  (EMI  =  -0.09, 
EM2  =  0.024,  EM3  =  0.026).  The  ebb  tidal  delta  is  shallow  and  symmetric  (EDI  = 
-1.7,  ED2  =  -0.34),  and  has  an  aera  of  1.02  square  miles. 

(6)  Cluster  6.  Cluster  6  consists  of  five  inlets,  again  geographi¬ 
cally  distributed  on  all  coasts.  These  inlets  are  the  widest,  longest,  and 
deepest  (W  =  4,100  feet,  I.  =  11,000  feet,  DMX  =  55  feet,  DMA  =  20  feet,  DCC  = 

26  feet,  EC1  =  -76,  EC2  =  42).  The  MIWC  is  shallow  and  only  moderately  asym¬ 
metric  ( EM 1  =  0.027,  EM2  =  0.028,  EM3  =  0.0003).  The  ebb  delta  is  relatively 
deep,  symmetric  and  large  (EDI  =  -1.2,  ED2  =  -0.79,  AED  =6.3  square  miles), 
though  not  as  deep  or  symmetric  as  some  other  inlet  clusters. 

(7)  Out  1 iers .  There  are  five  inlets  that  do  not  fit  any  of  the  clus¬ 
ters  defined.  Most  of  these  inlets  appear  relatively  large.  They  do  not  all, 
however,  cluster  together.  Hillsboro  appears  intermediate  to  a  combination  of 
clusters  1,  3,  and  4.  Doboy  and  Pass  Cavallo  are  close  to  each  other  as  are 
Willapa  and  Grays  Harbors.  Unfortunately,  for  inlets  of  this  size  it  is  often 
difficult  to  define  the  measures  such  as  crest  of  the  bar  and  length  of  channel 
because  of  either  the  complexity  and  size  of  the  inlet  or  poor  data. 

c.  Differences  Among  Clusters.  In  terms  of  width  and  length,  clusters  3 
and  4  are  the  most  similar,  with  cluster  5  inlets  slightly  wider  and  longer. 
Inlets  of  clusters  1,  2,  and  6  tend  to  be  much  wider  and  longer  than  clusters 
3,  4,  and  5.  The  primary  differences  between  clusters  3,  4,  and  5  are  in  depth 
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(Fig.  37, b)  and  relative  asymmetry  (EM2,  Table  5).  Cluster  3  tends  to  be 
shallower  than  cluster  4  which  is  shallower  than  cluster  5.  Clusters  1,  2, 
and  6  follow  a  similar  ordering;  however,  cluster  1  inlets  tend  to  be  some¬ 
what  shallower  than  cluster  4. 

Figure  37(e)  provides  an  interesting  summary  of  the  geometric  differences. 

In  this  figure  the  first  principal  components  of  inlet  cross-section  geometry 
(EMI)  (which  gives  relative  depth  of  the  cross  section  to  its  width)  are 
plotted  against  the  corresponding  component  of  ebb  delta  geometry  (EDI)  (which 
gives  relative  depth  of  the  delta  to  channel  depth).  The  arrangement  of  clus¬ 
ters  is  nearly  hierarchical.  Cluster  1  represents  the  shallowest  inlet  cross 
sections  (EMI  >  0)  and  has  correspondingly  a  relatively  deep  outer  bar  com¬ 
pared  to  channel  depths  (EDI  >  0).  Cluster  6  has  almost  as  shallow  an  MIWC, 
but  the  ebb  delta  is  much  shallower  compared  to  channel  depths  (EDI  <  0). 
Cluster  2  has  slightly  deeper  cross  section  than  cluster  1  or  6,  but  the  ebb 
delta  is  relatively  more  like  cluster  1.  Cluster  3  has  deeper  cross  sections 
than  clusters  1,  2,  or  6,  ar.d  a  shallower  bar  than  clusters  1  and  2.  Cluster 
5  has  a  much  deeper  cross  section  and  shallower  delta  than  those  clusters  pre¬ 
viously  mentioned.  Cluster  4  has  a  slightly  deeper  cross  section  and  shallower 
bar  than  cluster  5. 

If  the  clusters  are  ranked  in  ascending  value  of  mean  cross-sectional  areas, 
the  order  of  clusters  is  cluster  3  (2,000  square  feet),  cluster  4  (7,000  square 
feet),  cluster  5  (14,000  square  feet),  cluster  1  (29,000  square  feet),  cluster 
2  (40,000  square  feet),  and  cluster  6  (84,000  square  feet).  The  ranking  here 
follows  that  for  both  width  and  average  depth  as  should  be  expected.  It  again 
confirms  the  observation  that,  to  a  great  degree,  inlets  are  ordered  by  size. 

4.  Discuss i on. 

The  objective  of  the  analyses  just  presented  was  to  investigate  the  possi¬ 
bility  that  inlets  can  be  stratified  into  a  small  set  of  classes  on  the  basis 
of  inlet  geometry  and  to  describe  the  classes  found.  The  classification  anal¬ 
ysis  was  based  upon  the  multivariate  statistical  method  of  cluster  analysis. 

The  statistical  significance  of  the  clusters  found  was  tested  in  a  variety  of 
ways  and  showed  the  clustering  acceptable. 

The  analysis  determined  six  clusters  of  inlets  that  contained  all  but  five 
inlets  in  the  original  sample.  Inlets  in  the  individual  clusters  normally  in¬ 
clude  examples  from  two  or  three  coasts  which  must  imply  that  there  is  no  in¬ 
herent  reason,  based  on  the  geometric  parameters  analyzed,  for  stratification 
of  inlets  solely  on  a  geographic  basis.  The  inlet  clusters  can  be  arranged 
hierarchically  on  the  basis  of  size.  The  statistical  tests  give  strong  evi¬ 
dence  that  even  for  the  clusters  with  a  small  number  of  inlets  the  differences 
would  not  be  expected  by  sampling  errors  and  give  credence  to  a  hypothesis  of 
a  series  of  inlet  subpopulations. 

It  is  instructive  to  consider  further  the  question  of  whether  the  clusters 
represent  a  taxonometric  substructure.  The  discussion  on  t he  relationships 
between  the  geometric  parameters  provided  evidence  that  inlet  geometry  was 
organized  or  scaled  according  to  size.  The  clusters  are  likewise  organized. 

That  the  clustering  is  not  likely  due  to  a  lack  of  more  inlets  in  the  sample 
to  fill  in  the  gaps  can  be  exemplified  by  comparison  of  clusters  5  and  4.  The 
differences  in  mean  width  between  these  two  clusters  is  less  than  200  feet,  yet 
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the  me;in  DMA  for  cluster  3  is  almost  one-third  that  of  cluster  4.  This  would 
be  unexpected  if  the  differences  were  due  to  sampling  alone.  The  clusters  rep¬ 
resent  important  differences  in  inlet  geometry  other  than  just  a  scaling  process 
The  implications  of  this  in  relation  to  scaling  of  inlet  geometry  shown  previ¬ 
ously  are  considered  in  more  detail  in  Section  VII. 

The  clusters  presented  here  represent  the  taxonometric  structure  of  the 
sample  of  67  inlets  analyzed.  Because  it  is  felt  that  this  sample  is  fairly 
representative  of  the  range  of  inlets  on  the  U.S.  coasts,  the  statistics  de¬ 
rived  can  he  considered  estimates  of  the  population  statistics  for  all  inlets 
of  the  types  sampled  here.  As  has  been  indicated  by  the  outliers  in  the  analy¬ 
sis,  the  very  large  inlets  are  not  well  represented  in  this  analysis.  An  addi¬ 
tion  of  inlets  in  this  range  might  provide  additional  clustering.  Addition  of 
inlets  in  the  range  already  clustered  should  be  expected  to  redefine  the  sta¬ 
tistics  of  the  clusters  presented  but  should  not  force  combination  of  clusters 
already  defined. 

It  should  again  be  stressed  that  the  inlets  assigned  to  each  cluster  repre¬ 
sent  the  condition  of  the  inlet  at  a  particular  time  (Table  1),  not  necessarily 
today.  There  is  no  a  priori  reason  why  an  inlet  in  a  particular  cluster  must 
remain  there.  A  question  not  considered  here  is  the  stability  characteristics 
of  each  of  the  cluster  types. 

The  cluster  analysis  shows  that  inlet  geometry  has  a  strong  substructure 
rather  than  being  a  homogeneous  but  randomly  variable  population.  The  analysis 
is  based  on  a  sample  of  inlets  believed  typical  of  a  majority  of  U.S.  inlets 
and  clearly  implies  an  adjustment  of  inlet  geometry,  in  addition  to  the  scaling 
of  inlets,  according  to  size,  that  requires  explanation.  The  particular  clus¬ 
ters  presented  must  He  considered  a  first-order  stratification  of  inlets  on  the 
basis  of  inlet  geometry  that  could  be  refined  by  addition  of  more  inlets  or 
more  parameters.  The  small  number  of  clusters  provided  a  reasonable  framework 
for  investigating  the  geometric  variability  of  inlets.  It  is  intuitive  that  a 
finer  stratification  might  be  forced  througli  addition  of  more  parameters,  but 
it  is  perhaps  judicious  to  not  further  refine  the  classification  until  the 
physical  reasons  for  the  structure  presented  have  been  better  explained. 

V.  MATHEMATICAL  DEFINITION  OF  INLET  CLASSIFICATION 
1 .  Objectives. 

The  cluster  analysis  indicated  that  six  clusters  of  inlets  can  be  defined 
in  the  sample  analyzed.  If  the  view  is  taken  that  this  sample  provides  the 
estimates  of  the  multivariate  mean  and  variance  statistics  of  the  13  variables 
for  6  clusters  of  inlets,  it  is  desirable  to  form  a  series  of  equations  which 
mathematically  defines  the  classes.  What  is  sought  is  a  way  of  assigning  each 
member  of  the  sample  to  a  cluster  on  the  basis  of  a  probability  measure  and  a 
way  to  ultimately  derive  a  series  of  equations  which  allows  assignment  of  in¬ 
lets  not  in  the  analysis  to  the  clusters. 

The  method  used  to  produce  the  equations  is  a  discriminant  analysis.  The 
cluster  analysis  gave  an  indication  of  the  subpopulation  structure  of  inlet 
geometry  from  a  sample  of  inlets.  The  clusters  so  defined  were  used  to  develop 
estimates  of  the  subpopulation  mean  and  variance  statistics.  The  discriminant 
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analysis  takes  these  subpopulation  estimates  and  produces  equations  defining 
the  subpopulation  limits.  The  original  sample  is  reexamined  to  see  if  any 
misclassi ficat ion  has  occurred. 

2.  Mathematical  Considerations. 

Discriminant  analysis  is  a  widely  used  statistical  technique.  Any  number 
of  sources  are  available  for  reference  but  Kendall  and  Stuart  (1968),  Davis 
(1973),  and  Dixon  (1974)  are  particularly  helpful.  The  equations  presented 
here  are  from  Dixon  (1974).  They  are  presented  briefly  so  that  quantities 
later  presented  in  tables  can  be  more  readily  explained. 

The  basic  concept  behind  a  discriminant  analysis  can  be  seen  in  Figure  38 
if  it  is  assumed  that  two  clusters  (A,B)  have  been  found  in  analysis  on  two 
variables  Xj,  X2.  It  is  somewhat  intuitive  to  seek  a  line 

0  =  Q  +  AjXj  +  A2X2  (35) 

where  Q,  Aj  and  A2  are  constants,  in  (Xj,X2)  space  such  that  cluster  A  primarily 
lies  to  one  side,  with  B  on  the  other.  Depending  on  the  cluster  means  and  vari¬ 
ances,  there  may  be  many  or  few  elements  that  may  be  misclassified.  The  mis¬ 
classi  fication  occurs  either  because  random  errors  in  measurement  happen  or 
because  of  faulty  recognition;  i.e.,  an  element  of  A  is  incorrectly  called 
an  element  of  B.  The  line  defined  provides  a  basis  of  dividing  (Xj,X2)  space 
into  half  planes,  mathematically  defining  the  clusters  A  and  B.  The  equation 
of  the  line  is  defined  to  minimize  the  ratio  of  the  difference  between  the  pair 
of  multivariate  means  to  the  multivariate  variance  within  clusters. 


Figure  38.  Schematic  of  the  relationship  between  clusters 
A  and  B,  the  line  defined  by  the  discriminant 
analysis  for  variables  Xj  and  X2  by  which 
A  and  B  arc  classified. 
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The  problem  is  somewhat  more  complex  in  that  there  are  g(=  6)  clusters 
and  p(=  13)  variables.  The  simple  explanations  above  have  matrix  equation 
analogs.  The  equations  presented  here  are  defined  for  a  stepwise  discriminant 
analysis  where  variables  are  added  one  by  one  in  order  of  "best"  discrimination. 
This  is  to  imply  that  the  variable  added  into  the  analysis  is  that  variable 
which  maximizes  an  F  ratio  which  tests  the  hypothesis  that  the  multivariate 
mean  values  for  the  clusters  differ.  The  process  of  the  addition  of  variables 
continues  until  all  are  in  the  analysis.  At  any  step  in  the  analyses  consider 
that  1  <  r  <  p  variates  are  currently  in  the  analysis  and  that  p  -  r  are  not. 

The  following  parameters  are  needed.  Defining 


(a)  to  be  the  value  or  the 

inlet  from  the  m^’n  cluster 


variable  of  the 


(b)  Nw  number  of  inlets  in  cluster  m 

(c)  n  =  ni  +  n2  +  ...  +  n  :  total  number  of  inlets 

y 

_  \  8  nm 

(d)  X^  =  —  Z  Z  X^^:  the  mean  value  of  the  variable  X^ 

m=  l  k=  l 

_  i  nm 

(e)  rn^i  =  —  £  Xmki-  the  mean  of  X^  in  cluster  m 

nm  k=1 


,V/2 

:  the  standard 


smi  (n  -  1  ^  ^ mfei  ~ 

k=  i 


deviation  of  the  itn  variable  for  inlets  in  cluster  m 


^  ^  (*mki  ~  ~  m\ 

m=  l  k=l 

g  nm  —  _ 

(h)  (tfj)  =  £  £  i\nki  ~  i^mki  ~  %-j ) 

m= 1  k=  l 


the  matrices 


Wn  W12 

W21  W22 


T1  1  T 1 2 


T2 1  T22 


are  formed  where  Wn  and  Tu  are  r  *  r  cross-product  matrices 
(within  group  and  population-wide)  with  elements  w^j  and  t^j 

(g  and  h  above)  respectively,  for  the  r  variables  in  the  analysis 


at  the  given  step;  W12,  W21,  W22»  T12,  I21.  and  '*22  involve  p  -  r 
of  the  p  variates  not  currently  in  the  analysis.  (If  all  p  vari¬ 
ables  are  in  the  analysis,  W  =  Wjj  and  T  =  Tjj.) 

Forming  two  new  matrices 
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g  discriminant  functions  can  be  formed,  one  for  each  cluster. 
The  g  discriminant  functions  are  defined  by 


(44) 


(45) 


r 

Sj imk  =  C mo  +  ^  C mj  (46) 

j  =  l 

where  is  the  value  for  the  inlet  of  the  cluster  on  the  m 

(of  g)  discriminant  functions.  The  coefficients  in  the  equation  are  defined  as 

h 

C mi  =  (n  ~  g)  aij  (4  7) 

j=l 

and 

1  r 

^mo  ~  2  ^  CffU  (48) 

j  =  l 

where  r  is  the  number  of  variables  in  the  analysis  at  the  given  step  and  the 
other  variables  are  as  previously  defined. 

Given  two  arbitrary  clusters  (m,£)  the  Mahalanobis  squared  distance,  D2 
is 


r 

L 

i=  1 


(C *ri  -  Cj ii) 


{m^i  ~  £.Xf) 


(49) 


This  function  measures  the  statistical  difference  between  the  multivariate 
means  for  the  two  clusters.  An  F  value  for  the  difference  between  clusters 
for  testing  purposes  is 


F  =  (n  -  g  -  r  +  1) 
ml  "  r(n  -  g)  (n^  +  n-.) 


<1 


(50) 
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with  r  and  n  -  g  -  r  +  1  degrees  of  freedom.  If  the  value  of  F  is  signifi¬ 
cant  at  some  level  a,  then  the  chance  that  the  multivariable  means  for  the 
two  clusters  are  equal  is  no  more  than  a.  The  values  of  a  normally  used  are 
5  and  1  percent.  If  the  F  test  is  significant  at  one  of  these  a,  the  clusters 
are  considered  distinct. 

It  is  also  possible  to  test  the  equality  of  all  group  means  simultaneously 
by  forming  the  ratio 

detfWn) 

u  =  - -  (5  1 J 

det (Til) 

and  calculating 

F  =  1  ~  u1/2  x  yz  +  1  -  r(g  -  l)/2  (52) 

u1^2  r(g  -  1) 

where  z  =  r2(g  -  l)2  -  4/r2  +  (g  -  l)2  -  5  if  r2  +  (g  -  l)2  ^  5,  z  =  1  other¬ 
wise,  and 

y  •  n  -  (53) 

This  F  ratio  has  r(g  -  1)  and  yz  +  1  -  r(g  -  l)/2  degrees  of  freedom.  The  F 
test  again  tests  whether  or  not  the  clusters  defined  are  considered  distinct 
among  each  other  on  some  level  of  significance. 

Finally,  it  is  possible  to  calculate  the  posterior  probability  that  inlet 
k  of  cluster  l  actually  belonged  to  cluster  m; 

Pra  exp(SJDr?P 

P^  =  - - (54) 

Z  Pz  exP 
i=  1 

with  defined  previously  and 

nn 

P,  .  /  (55) 

the  prior  probability  that  the  inlet  belonged  to  a  given  cluster,  q.  As  dis¬ 
cussed  in  the  discriminant  analysis  of  the  data,  equation  (54)  can  be  used  to 
calculate  the  probability  that  an  inlet  not  in  the  analysis  belongs  to  a  par¬ 
ticular  cluster. 

3.  Discriminant  Analysis  of  Inlet  Data. 


A  stepwise  discriminant  analysis  as  described  was  performed  incorporating 
all  13  variables.  Then  a  more  restricted  analysis  was  performed  involving  only 
the  six  variables  DMX,  DMA,  W,  DCC,  L,  and  AliD  to  see  if  simpler  discriminant 
functions  could  be  derived. 


a.  Stepwise  Analyses.  In  the  first  analysis  the  objective  is  to  find  what 
minimum  set  of  variables  provided  adquate  cluster  discrimination;  i.e.,  what 
set  of  variates  produces  a  series  of  equations  that  separates  the  majority  of 
the  inlets  into  the  original  clusters  derived  on  the  basis  of  equation  (49)  and 
listed  by  equation  (50).  For  this  to  occur,  the  F  ratio  must  be  significantly 
above  a  given  level,  taken  to  be  1  percent  in  this  study,  for  every  pair  of  the 
original  six  clusters.  The  analysis  proceeded  by  selecting  one  by  one  from  the 
set  of  13  variables  the  ones  that  increase  discrimination.  Using  a  1-percent 
level  of  significance,  the  first  three  variables  added  were  EM2,  DCC,  EM3,  and 
provided  acceptable  discrimination  at  1  percent.  The  discriminant  functions 
derived  are  given  in  Table  6  as  are  the  F  ratios  and  a  classification  matrix 
which  indicates  how  many  inlets  of  each  cluster  are  classified  as  belonging  to 
another  cluster.  The  coefficients  of  the  discriminant  function  given  in  Table 
6(a)  can  be  used  to  calculate  the  probabilities  that  inlet  belongs  to  a  given 
cluster  using  equation  (54).  The  matrix  of  F  ratios  (Table  6,b)  is  used  to 
test  the  significance  that  any  two  clusters  are  statistically  different  based 
on  a  critical  value  of  the  F  ratio.  In  this  case,  the  F  ratio  must  have  a 
value  of  at  least  4.3  in  order  to  accept  the  hypothesis  that  the  multivariate 
means  for  pairs  of  clusters  are  not  equal.  The  classification  matrix  (Table 
6,c)  indicates  the  number  of  inlets  in  each  cluster  classified  as  being  in 
given  clusters  based  on  the  discriminant  functions.  The  classification  matrix 
indicates  that  a  total  of  11  of  the  62  inlets  analyzed  (the  other  five  were 
neglected)  are  misclassif ied;  i.e.,  on  the  basis  of  the  functions  developed, 
the  inlets  are  in  the  wrong  cluster. 


Table  6.  Descriminant  analysis  results  for  three  variables  (DCC,  EM2, and  EM3) 

1  a.  Coefficients  for  discriminant  functions  based  on  three  variables. _ 
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Classification  matrix. 


Number  of  cases  classified  into  group 


Croup  Type  1  Type  2  Type  3  Type  4  Type  5  Type  6 


3  0  0  0 

3  3  0 

2  0  2 

0  0  0 

4  0  0  0  8 

8  1  0  0  0 

3  and  51  degrees  of  fr'cdom  and  an  1  value  of  4.3  is  significant  at  a  l-perccnt  level. 


The  next  step  was  to  continue  adding  variables  in  order  to  determine  what 
minimal  group  of  variables  provides  adequate  cluster  discrimination  and  is 
closest  to  the  initial  clusters  determined.  After  nine  variables  (EM2,  DCC, 
DM3,  DMA,  EMI,  HD2,  EDI,  EC2,  and  W)  were  added,  the  cluster  discrimination 
did  not  differ  in  classification  matrix  for  an  analysis  with  the  full  set  of 


variables.  In  this  instance  only  one  inlet  is  misclassi fied.  This  is  Tilla¬ 
mook  of  cluster  6  which  is  classified  in  cluster  2  by  the  analysis.  The  dis¬ 
criminant  functions  for  the  full  case  of  13  variables,  the  F  ratios  for  cluster 
pairs,  and  the  classification  matrix  are  given  in  Table  7.  Table  8  gives  the 
posterior  probabilities  for  each  of  the  62  inlets;  i.e.,  the  probability  based 
on  the  analysis  that  the  inlet  should  be  in  a  given  cluster.  Again,  the  sta¬ 
tistical  tests  on  the  cluster  pairs  are  significantly  above  1  percent. 

b.  Truncated  Parameter  Analysis.  The  second  discriminant  analysis  was 
performed  only  with  the  variables  DMX,  DMA,  W,  DCC,  L,  and  AED.  It  was  rec¬ 
ognized  that  these  variables  are  more  easily  measured  than  the  other  variables. 
Thus,  if  these  variables  provide  a  viable  cluster  discrimination,  they  would 
provide  a  simpler  basis  for  classification. 


Table  9  provides  the  F  ratios  for  cluster  pairs  and  the  classification 
matrix  for  this  analysis.  For  the  F  ratios  computed  with  \>j  =  6  and  v2  =  51 
degrees  of  freedom,  a  value  greater  than  3.25  is  required  to  reject  the  hypoth¬ 
esis  of  equal  multivariate  means  (or  to  accept  the  hypothesis  that  the  clusters 
are  different).  Table  9  indicates  that  numerous  cluster  pairs  (1-3,  1-4,  2-4, 
2-5,  3-4,  4-5)  do  not  meet  this  criteria.  Hence,  the  simplification  to  these 
few  variables  does  not  reproduce  the  clusters  originally  formed  on  the  basis  of 
all  13  variables. 


Table  7.  Discriminant  analysis  results  for  13  variables. 
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The  V  test  has  13  ind  44  degrees  of  freedom  and  an  F  value  of  2.6  is  significant  at  a  1-pcrcent  level 
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Table-9.  Descriminant  analysis  results  for  six  variables  (DMX,  DMA,  W,  DCC,  L,  and  A£l)) . 
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a . 

Coefficients  for  discriminant  functions 

based  on  six  variables. 
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Classification  matrix. 

Number  of  c..ses  classified  into  group 
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4.  Discussion. 


Since  it  is  not  possible  to  discriminate  the  original  clusters  with  a  trun¬ 
cated  parameter  set,  the  discriminant  functions  must  be  chosen  from  the  step¬ 
wise  analysis  of  the  complete  parameter  set.  It  would  be  possible  to  use  the 
three  variable-based  functions  (Table  6).  However,  these  variables  are  rela¬ 
tively  difficult  to  derive.  It  is  recommended  that  the  functions  based  on  the 
full  set  (Table  8)  be  used.  The  effort  involved  is  not  that  much  greater  than 
for  the  th^ee  variable  cases. 

Several  implications  result  from  the  discriminant  analysis.  First,  the 
stepwise  discriminant  analysis  again  offers  strong  evidence,  based  upon  multi¬ 
variate  analyses  of  variance,  that  the  clusters  are  well  defined.  The  misclas- 
sification  of  only  one  inlet  and  the  strength  of  the  posterior  possibilities 
(Table  9)  provide  "evidence  of  this.  The  discriminant  functions  derived  provide 
a  basis  of  assignment  of  the  inlets  to  the  original  clusters. 

The  second  implication  seen  is  the  dependence  of  the  analysis  on  the  shape 
functions  introduced  (liMl,  LM2,  FM3,  FC2,  liDl,  and  HD2) .  The  clusters  cannot 
be  discriminated  without  them.  This,  in  conjunction  with  their  appealing  geo¬ 
metric  interpretation,  suggests  that  they  perform  an  adequate  and  useful  repre¬ 
sentation  of  inlet  geometry. 

Thirdly,  it  is  seen  that  DMX,  AliD,  LC1,  and  L  are  perhaps  redundant  param¬ 
eters  because  in  the  stepwise  discrimination  they  provide  no  further  refinement 
They  are  maintained  for  purely  descriptive  purposes. 


through  the  discriminant  analysis,  functions  have  been  derived  that  mathe¬ 
matically  define  the  clusters.  From  these  functions  it  is  possible  to  derive 


the  probability  that  an  inlet  belongs  to  a  given  cluster.  The  relationships 
and  cluster  discrimination  appear  statistically  significant  at  a  very  high 
leve  1 . 


VI.  EVALUATION  OF-'  THE  PROBABILITY  THAT  AN  INLET 
NOT  IN  THE  ANALYSIS  BELONGS  TO  A  CLUSTER 

The  analyses  presented  in  this  report  represent  a  classification  of  67  in¬ 
lets  at  selected  times.  It  is  desirable  to  have  a  method  for  placing  an  inlet 
not  in  the  analysis  into  the  classification.  Such  a  method  would  also  allow 
comparison  of  the  later  condition  of  an  inlet  in  the  analysis  with  the  condi¬ 
tion  originally  used.  In  the  strictest  sense,  the  appropriate  method  would  be 
to  enter  the  new  inlet  into  the  analysis  and  redo  the  entire  classification  and 
discriminant  analysis.  This  is  a  time-consuming  and  laborious  exercise.  It 
should  not  be  considered  unless  a  large  number  of  inlets  arc  to  he  entered. 

The  following  method  represents  a  simplified  method  for  estimating  the  proba¬ 
bility  that  the  new  inlet  belongs  to  a  given  cluster.  The  principal  assump¬ 
tions  required  are  that  f a J  the  new  inlet  does  belong  to  one  of  the  clusters 
and  (b)  the  discriminant  functions  calculated  in  Section  V  are  adequate  sta¬ 
tistical  formulas  describing  the  classification. 

1 .  Data  Preparation. 

To  use  the  discriminant  functions,  it  is  necessary  to  measure  and  calculate 
the  13  variables  used  in  the  same  fashion  as  they  were  for  this  report.  for 
the  variables,  DMX,  DMA,  IV,  DCC,  L,  and  Aid),  the  definitions  are  straightfor¬ 
ward.  l-'or  the  variables  EMI,  EM2,  EM3,  EC1,  EC2 ,  EDI,  and  EI)2,  the  definitions 
are  less  straightforward.  The  following  discussion  describes  how  the  values 
can  be  calculated  efficiently. 

Taking  the  minimum  inlet  width  cross-section  eigenvectors  (e_j ,  e2,  e_3)  as 
an  example,  the  following  steps  from  Section  II  are  necessary  to  find  the 
weightings  EMI,  EM2,  and  EM3: 

(a)  A  cross  section  at  the  minimum  inlet  width  is  drawn. 

(b J  The  60  evenly  spaced  points  across  the  profile  are  located 
and  the  depth  recorded  with  the  order  in  which  the  depths  are  listed 
based  on  the  convention  of  Section  II.  Each  depth  is  divided  by  W. 

(c)  The  mean  depth  and  the  standard  deviation  for  each  depth 
location  (calculated  and  listed  in  App.  D)  are  used  to  normalize 
each  depth 


where  dj  is  the  newly  measured  depth,  d^  and  S^  are  the 
mean  and  standard  deviation  respectively,  and  d;-  is  the 
normalized  value. 


(d)  A  new  vector  D  =  (d*,  d*  ...,  d*  )  is  formed 
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(e)  For  ei  the  coefficients  are  taken  from  Appendix  E  to  form 
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ing  the  procedures  of  Section  II  and  using  values  of  Appendixes  D  and 
E.  A  final  check  would  be  to  find  an  inlet  similar  in  characteristics 
to  one  in  this  report  and  compare  values  to  see  if  the  ones  newly  cal¬ 
culated  look  reasonable. 


2 .  Probability  Calculations. 

To  calculate  the  desired  probability,  it  is  necessary  to  evaluate  the 
following  values,  where  i  will  range  from  1  to  6. 

13 

%  =  C i0  +  L  C ij  Xj  (62) 

j  =  l 

where  C'£0  is  the  constant  on  the  i ^  discriminant  function  (eq.  48),  C^j 
is  the  coefficient  for  the  j ^  of  13  variables  for  the  i ^  discriminant 
function  (eq.  47),  and  X*  is  the  value  of  the  j  ^  variable.  The  subscript 
i  relates  the  discriminant  function  to  the  appropriate  inlet  group  and  equa¬ 
tion  (62)  is  seen  to  be  analogous  to  equation  (46).  In  actual  computations, 
the  values  of  C-  0  and  C7-j  would  come  from  Table  7(a)  . 

The  computation  of  the  six  values  gives  only  the  values  of  the  dis¬ 

criminant  functions.  The  probability  that  the  inlet  belongs  to  inlet  group  i 
is  p,;  calculated  by 


*  =  vi  cxp  (s?;) 

i  6 

7.  p^  exp  (S.) 

J  ° 


(63) 
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where 


P  • 

J 


is  the  a  priori  probability  of  belonging  to  a  given  inlet  group 


(64) 


where  Nr  is  the  number  of  inlets  in  the  j th  group  (N  ■  =  31,  6,  4,  4,  12,  5 

d 

for  j  =  1  to  6,  respectively).  Liquation  (63)  is  analogous  to  equation  (54). 


These  equations  for  calculating  the  discriminant  functions  and  the  proba¬ 
bilities  iiave  been  programed  into  a  simple  computer  program  given  in  Appendix 
F.  It  is  written  in  a  time-sharing  format  with  unformatted  input  and  output. 
With  minor  modifications,  it  can  lie  operational  on  almost  any  digital  computer. 

3.  A  Simple  i.xamp  le . 

Fable  10  provides  a  simple  example  of  input  to  the  program-the  values  of 
tLie  13  variables  for  a  hypothetical  inlet.  The  results  indicate  that  the  inlet 
group  for  which  the  probability  is  highest  is  group  3  with  a  value  of  0.85.  The 
second  highest  is  group  4  with  a  value  of  0.14. 


Fable  10.  A  simple  example  of  the  calculation  of  the  proba¬ 
bility  that  an  inlet  belongs  to  a  given  group. 


Input  to  computer  program  (App.  F) 


DUX  -25. 0 

UCC  -4.0 

EMI  -0.1 
|  EC1  20.0 

I  EDI  -3.0 

DMA  -15.0 

L  1,500.0 
*EM2  0.08 

EC2  -1.0 

ED2  1.0 

N  1,000.0 

EM3  0.0 

AED  0.15 

1  Output  from  program 

Group 

Si 

Pi 

Pi 

(a 

priori  probability) 

(probability  from 

discriminant  analysis) 

1 

34.3 

0.500 

0.000 

2 

23.5 

0.097 

0.000 

3 

43.0 

0.065 

0.853 

4 

46.2 

0.065 

0.  141 

5 

41.8 

0.  193 

0.001 

6 

-2.9 

0.081 

0.000 

1EM2  input  as  absolute  value. 


VII.  IN  LIFT  GEOMETRY:  A  SUMMARY  DISCUSSION 

The  initial  objective  of  the  study  was  the  definition  of  a  series  of  param¬ 
eters  that  satisfactorily  describes  the  major  components  of  the  inlet  throat  and 
outer  bar.  After  preliminary  studies  13  variables  were  selected  incorporating 
its  factors  that  describe  both  the  physical  dimensions  of  the  inlet  and  the  shape 
(or  geometry).  Many  of  the  parameters  have  been  used  previously  and  have  anal¬ 
ogous  roles  in  describing  the  hydraulic  character  of  inlets.  Other  parameters 
are  new  and  arc  a  result  of  recent  research  into  the  field  of  shape  analysis. 

The  results  of  the  statistical  analyses  performed  in  the  study  suggest  that 
the  parameters  chosen  are  a  set  sufficient  for  resolving  a  series  of  principal 
questions  regarding  the  variability  of  inlet  geometry.  Through  careful  measure¬ 
ment  of  these  variates  a  series  of  relationships  among  the  geometric  parameters 
has  been  defined. 
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The  eigenvector-generated  description  of  geometric  components  of  inlet  bar 
and  channel  shapes  provided  a  clear,  concise  indication  of  the  covariability  of 
these  factors  with  the  other  variates.  Although  the  mathematics  of  their  deri¬ 
vation  is  not  trivial,  the  resulting  parameters  can  be  interpreted  in  a  simple 
fashion.  As  the  discriminant  analysis  indicates,  they  have  an  important  role 
in  defining  the  clustering  found.  It  is  expected  that  these  parameters  and 
others  like  them  will  find  an  increasingly  important  role  in  the  description 
and  analysis  of  land forms. 

The  set  of  parameters  chosen  and  the  particular  definitions  used  are  by  no 
means  the  only  ways  of  representing  inlet  morphology.  It  is  evident  that  some 
of  the  parameters  are  redundant,  but  these  have  been  retained  for  descriptive 
purposes.  Any  future  classification  analyses  should  consider  their  elimination. 
One  major  variate  not  directly  used  in  the  study  is  ebb  delta  volume  defined 
along  the  lines  suggested  in  Dean  and  Walton  (1975J.  It  is  obvious  that  some 
of  this  information  is  retained  in  the  ebb  delta  eigenvectors;  the  principal 
difficulty  lies  in  the  consistent  definition  of  the  base  surface  above  which 
the  volume  is  taken.  In  areas  with  multiple,  and  possibly  overlapping  inlet 
deltas  as  is  the  case  in  a  number  of  coastal  areas,  the  base  surface  is  diffi¬ 
cult  if  not  impossible  to  define. 

One  principal  goal  of  the  classification  study  was  to  analyze  the  relation¬ 
ships  among  the  geometric  parameters  chosen.  It  was  a  reasonable  expectation 
to  find  moderate  partial  correlations  among  the  parameters.  A  major  outcome 
of  the  study  is  the  scries  of  strong  relationships  found  that  is  statistically 
significant  at  high  levels.  It  is  apparent  from  the  study  that  many  of  the 
geometric  factors  are  interrelated  in  a  predictable  manner  and  that  the  param¬ 
eters  DMX,  DMA,  DCC,  I,,  and  /Mi D  in  particular  can  he  co-rclatcd.  Likewise, 
when  inlet  width  did  not  appear  as  a  controlling  parameter  in  the  variable 
cross  plots,  the  strength  of  the  relationships  of  so  many  parameters  to  the 
area  of  the  minimum  inlet  width  cross  section,  A^  (estimated  by  W  *  DMA), 
was  unexpected.  These  relationships  to  A,,  provided  a  key  to  an  understanding 
of  the  adjustment  of  inlet  geometry. 

The  major  i  mp  1  i  cat  i  on  of  t  lie  relationships  found  is  the  large  scale  coadjust¬ 
ment  of  inlet  geometry  that  appears  scaled  by  the  parameter  Ac.  As  Ac.  in¬ 
creases,  the  inlet  channel  becomes  less  incised  into  the  outer  bar,  the  relative 
depth  of  the  channel  across  the  bar  increases,  the  channel  length  is  increased, 
the  ebb  delta  area  enlarges,  and  channel  depths  deepen.  It  is  apparent  that 
this  adjustment  occurs  as  a  response  to  the  wave  and  tide  processes  and  that 
the  relationships  found  are  statistical  summaries  of  complex  hydraulic-sediment 
interactions  in  widely  ranging  geologic  settings.  An  important  component  of 
future  research  should  be  an  effort  to  place  the  relationships  in  the  perspec¬ 
tive  of  the  wave  and  tide  processes  at  the  inlets  studied. 

The  relationships  found  may  have  direct  design  implications  for  a  number  of 
practical  engineering  problems  involving  the  design  of  inlet  modifications.  If 
the  tidal  prism-minimum  cross-sectional  area  relationship  of  O'Brien  (19.fl)  is 
assumed,  it  is  apparent  that  given  one  prism,  then  one  A,.,  and  hence  a  narrow 
range  of  inlet  geometries  is  possible.  Whether  this  is  in  fact  realistic,  or 
whether  the  geometry  predicted  is  an  equilibrium-type  form  that  the  condition 
might  produce  is  unknown.  A  further  implication  of  the  analysis  is  that  the 


internal  adjustment  of  DMX,  DMA,  1,,  and  DCC  (among  other  parameters)  suggests 
that  the  detail  of  the  geometry  is  highly  coordinated,  liven  the  shape  of  the 
channel  profile  appears  to  he  determined  to  a  large  degree.  The  adjustments 
established  in  this  report  result  in  inlet  channel  geometries  unlike  those 
currently  used  in  design  practice,  lixamples  of  this  construction  of  constant 
depth  channels  while  most  of  the  natural  channels  have  a  significant  slope  and 
construction  of  steep-sided  or  U-shaped  channel  cross  sections  while  natural 
channels  are  more  gently  sloping.  ITie  implication  of  this  upon  the  success  or 
failure  of  the  proposed  modifications  to  an  inlet  is  unknown. 

It  is  not  clear,  if  the  relationships  developed  here  are  applied  in  the 
design  of  a  jettied  or  dredged  channel,  how  they  are  to  be  applied,  or  whether 
they  will  be  a  valid  prediction  of  inlet  response.  As  an  example,  if  an  inlet 
is  modified  to  have  a  highly  constricted  cross-sectional  area,  it  is  not  evi¬ 
dent  that  DMX,  DCC,  or  L  must  respond  as  predicted  because  the  throat  geometry 
inherently  associated  with  the  given  cross  section  is  not  necessarily  preserved 
in  the  proposed  modification.  it  is  perhaps  reasonable  to  expect  that  the  equa¬ 
tions  will  produce  bounds  for  the  response;  however,  even  this  is  not  assured. 

It  is  thus  recommended  that,  before  these  relationships  are  used  in  design,  more 
research  and  experimental  effort  be  given  to  provide  a  better  understanding  of 
their  implications. 

The  second  principal  goal  of  the  study  was  to  investigate  a  possible  classi¬ 
fication  of  inlets  based  upon  inlet  geometry.  The  objectives  were  to  see  if  a 
subpopulation  structure  existed  and,  if  it  existed,  to  define  it.  The  result¬ 
ing  classification  would  produce  a  better  understanding  of  inlet  variability 
and  would  aid  personnel  in  design  projects  to  find  prototype  inlets  of  similar 
characteristics.  The  results  of  both  the  cluster  analysis  and  the  discriminant 
analysis  indicate  the  presence  of  at  least  six  well-defined  clusters  or  types 
of  inlets  based  on  geometry.  The  cluster  analysis  provided  the  taxonometric 
structure  of  the  set  of  inlets  analyzed  and  produced  a  preliminary  classifica¬ 
tion.  The  discriminant  analysis  further  refined  the  classification  to  produce 
a  series  of  functions  that  allows  assignment  of  inlets  to  the  six  clusters.  It 
was  also  recognized  that  the  very  large  inlets  were  not  represented  sufficiently 
in  the  analysis  to  provide  reliable  clustering  or  discrimination. 

hxamination  of  the  relationships  among  the  clusters  derived  indicated  that 
the  scaling  process  evident  in  the  other  analyses  is  preserved.  The  signifi¬ 
cant  result  of  the  classification  analyses,  however,  is  that  the  clusters  sys¬ 
tematically  organize  inlets  on  the  basis  of  width,  depth,  and  shape  in  a  way 
that  accounts  for  some  of  the  scatter  observed  in  simple  relationships  among 
parameters.  Thus,  clusters  3  and  4  differ  only  slightly  in  width,  but  greatly 
in  terms  of  the  other  variables.  Clusters  4  and  5  have  similar  depths  and 
widths  of  the  cross  section  but  differing  depths  at  the  crest  of  the  outer 
bar,  and  so  forth  for  other  pairs  of  clusters.  In  part,  the  clusters  would 
appear  to  account  for  some  of  the  variability  in  geometry  that  may  be  attrib¬ 
utable  to  wave  action.  Whether  this  is  indeed  the  case  will  require  future 
research  relating  the  wave  and  tide  processes  to  the  geometry. 

To  summarize  the  implications  of  both  the  classification  and  parameter 
variation  studies,  it  is  essential  to  realize  that  they  are  complementary, 
together  explaining  in  some  detail  the  systematic  organization  of  tidal  inlet 
geometry.  The  cross-sectional  area,  A,, ,  can  be  simpl i st ical ly  considered  as 
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a  scaling  parameter  that  relates  the  relative  magnitudes  of  wave-  and  tide¬ 
generated  sediment  transport.  With  increasing  area,  the  tide  processes  appear 
to  dominate  the  wave  processes;  inlet  geometry  varies  accordingly.  The  size  of 
the  cross-sectional  area  scales  the  bar  and  channel  geometry  in  a  fairly  reg¬ 
ular  pattern.  As  a  result  it  must  be  concluded  that,  except  for  small  inlets, 
the  t i de-gene  rated  processes  to  a  large  degree  determine  the  geometry,  liven 
in  the  smaller  inlets  tidal  control  is  still  evident. 

The  classification  analyses  substantiate  the  scaling  relationship,  but 
indicate  that  there  are  systematic  deviations  of  inlet  geometry  not  fully 
explained  by  the  scaling  relationship.  The  clusters  found  represent  this 
organization.  It  is  clear  also  that  if  the  deviations  away  from  the  scaling 
relationships  were  purely  random  it  would  be  fortuitous  to  have  the  discrimi¬ 
nant  analysis  significant  at  a  high  level.  The  implication  must  be  that  the 
organization  is  real  and  forced  by  some  underlying  cause. 

The  absence  of  width,  and  less  so  average  depth,  as  important  scaling  fac¬ 
tors  in  many  relationships  shown  in  Section  III  and  the  absence  of  relation¬ 
ships  between  LC2  and  EM2  and  other  parameters  underscore  the  lack  of  a 
parameter  that  performs  a  scaling  for  wave  action  in  the  way  that  ka  scales 
the  larger  scale  geometry.  hC2,  which  describes  the  slope  of  the  channel 
profile,  is  probably  related  to  onshore-offshore  sediment  transport  by  waves 
at  the  edges  of  the  delta  and  perhaps  less  so  to  longshore  transport.  liM2 , 
which  gives  cross-section  asymmetry,  would  be  more  likely  related  to  wave-caused 
longshore  transport.  If  this  is  so,  there  would  not  be  necessarily  a  good  rela¬ 
tionship  between  IiM2  and  LiC2 . 

In  addition  to  meeting  the  primary  objectives  of  the  study,  the  parameters 
chosen  and  the  analyses  performed  place  the  variation  of  inlet  geometry  into 
a  better  perspective  of  the  varying  influences  of  waves  and  tides.  It  is  im¬ 
portant  to  recognize  that  the  conclusions  drawn  are  based  upon  an  interpreta¬ 
tion  of  the  morphology  of  tidal  inlet  systems  and  by  design  have  not  involved 
estimates  of  the  hydrodynamic  processes.  This  approach  was,  in  part,  followed 
to  see  if  a  natural  organization  of  inlet  variability  was  evident  which  would 
motivate  research  into  the  correspondence  between  process  and  form.  The  results 
of  the  study  justify  the  need  for  performance  of  this  work. 

As  is  often  the  case  in  studies  of  this  type,  more  questions  are  generated 
than  resolved.  In  particular  the  results  should  motivate  study  into  the  rela¬ 
tionships  between  the  types  of  inlet  geometry  found  and  tidal  prisms,  ranges 
and  currents,  net  and  gross  longshore  drift,  and  onshore  and  offshore  sediment 
movement  by  waves.  The  certainty  of  whether  the  inlet  types  are  natural  by¬ 
passers  of  sand  or  not  needs  to  be  established.  The  variation  of  inlet  sta¬ 
bility  by  inlet  type  requires  examination.  Finally,  the  number  of  inlets 
analyzed  should  be  increased. 

Several  of  the  parameters  and  a  number  of  analyses  used  in  this  report  have 
not  been  extensively  used  in  cither  the  engineering  or  geologic  literature. 

When  the  scientific  questions  of  a  study  involve  multidimensional  variation, 
there  is  a  wealth  of  statistical  procedures  that  can  be  used  in  a  rigorous 
method  of  investigation.  Given  the  complex  multidimensional  variation  typical 
of  inlet  geometry,  it  is  difficult  to  sec  how  the  results  obtained  here  could 
have  been  achieved  using  only  one-dimensional  methods. 
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VIII.  SUMMARY 


Parameters  have  been  devised  that  measure  and  describe  inlet  throat  and 
ebb  delta  morphology.  When  these  parameters  are  measured  for  inlets  where 
sufficient  chart  data  are  available,  the  parameters  are  shown  to  vary  in  a 
consistent  fashion  that  appears  to  be  scaled  according  to  the  relative  magni¬ 
tude  of  the  tidal  processes.  Wien  tiie  subpopulation  structure  is  examined, 
inlets  can  be  initially  clustered  into  six  classes  which  can  be  mathematically 
discriminated.  The  classification  provides  a  systematic  organization  of  inlet 
geometry  that  is  related  to  deviations  from  the  basic  scaling  relationship 
probably  due  to  the  influence  of  wave  action.  The  relationships  and  classi¬ 
fication  found  are  statistically  significant  at  high  levels  which  provides 
confidence  in  the  results. 

In  a  statement  generally  attributed  to  O'Brien,  there  are  said  to  be  two 
types  of  inlets:  large  and  small.  This  study  to  a  large  degree  confirms  this 
observation,  but  beyond  that  shows  the  adjustment  of  inlet  geometry  to  be 
systematic  and  predictable. 
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APPENDIX  D 


MEAN  AND  STANDARD  DEVIATION  VECTORS  REQUIRED  IN  EIGENVECTOR  ANALYSIS 


a.  Normalized  cross  section 
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APPENDIX  E 

EIGENVECTORS  FOR  THE  MINIMUM  INLET  CROSS  SECTION, 
CHANNEL  PROFILE,  AND  EBB  DELTA  GEOMETRY  ANALYSES 

a.  Minimum  inlet  cross  section 


EMI 


EM2 


EM3 


b.  Channel  profile 


c.  Ebb  delta  geometry. 
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54 

0.1781 
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APPENDIX  F 

COMPUTER  PROGRAM  FOR  THE  CALCULATION  OF  THE  PROBABILITY 
THAT  AN  INLET  BELONGS  TO  AN  INLET  GROUP 

This  appendix  contains  a  simple  computer  program  that  incorporates  the  dis¬ 
criminant  functions  of  Table  7(a)  and  can  be  used  to  calculate  the  posterior 
probability  that  an  inlet  belongs  to  a  group.  The  program  is  written  in  time¬ 
sharing  FORTRAN  for  a  Honeywell  637  computer.  Input  and  output  are  free  field. 
A  person  with  modest  programing  capabilities  should  be  able  to  make  the  program 
run  on  any  FORTRAN  system. 

The  following  input  to  the  program  is  for  one  inlet: 

(a)  DMX,  DMA,  W 

(b)  DCC,  L 

(c)  EMI,  EM2,  EM 5 

(d)  EC1,  EC2 

(e)  EDI,  ED2,  AED 

EM2  must  be  input  as  an  absolute  value. 

The  output  is  the  values  of  the  discriminant  function,  the  prior  proba¬ 
bility,  and  the  posterior  probability  as  a  triplet  for  each  inlet  group. 
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100  CONTINUE 

CRLCULHTE  THE  F'PDFRF  IL I  TV  FRIED  ON  R  PRIORI  RND 
DI scpiminhnt  function: 

:  i  ::=  o. 

DM  r.5  11=1,5 
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Dll  )t.h  II  =1«5 

i , . ii  p:trp'  i*  ••-p< u  "*e::P'  ::  ■  it  o/sf: 

I'll  .300  1  =  1,5 
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